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Abstract 

The Hamilton Jacobi formalism generalized to 2 dimensional field theories ac- 
cording to Lepage's canonical framework is applied to several relativistic real 
scalar fields, e.g. massless and massive Klein Gordon, Sinh &: Sine Gordon, 
Liouville and 4> 4 theories. The relations between the Euler Lagrange and the 
Hamilton Jacobi equations are discussed in DeDonder and Weyl's and the 
corresponding wave fronts are calculated in Caratheodory's formulation. 
Unlike mechanics we have to impose certain integrability conditions on the 
velocity fields to guarantee the transversality relations and especially the dy- 
namical equivalence between Hamilton Jacobi wave fronts and families of ex- 
tremals embedded therein. 

Backlund Transformations play a crucial role in solving the resulting system 
of coupled nonlinear PDEs. 



1 Introduction 



Varying a relativistically invariant action integral leads to covariant Euler-Lagrange 
equations. However, if one wants to reformulate the theory in terms of the conven- 
tional canonical Hamiltonian framework one has to break the manifest covariance by 
distinguishing a time variable and regarding the other "spatial" coordinates as "indices" 
representing an infinite number of degrees of freedom. The method is widely known from 
elementary particle physics, canonical gravity and other field theories. This approach, 
however, can obscure a part of the rich geometrical structure contained in a generally 
covariant framework, at least on the classical level ffl: 
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Utilizing Cartan's theory of alternating forms Lepage and others showed that a large va- 
riety of algebraically inequivalent covariant Hamiltonian formulations, including Hamil- 
ton-Jacobi equations, exists, e.g. that of DeDonder & Weyl ]1U) (Tl[ and that of Cara- 



theodory |12|], where only the latter one provides a Hamilton-Jacobi equation the as- 



sociated wave fronts of which have the same nice transversality properties with respect 
to the extremals as one has in mechanics. Reformulations and special examples of this 
covariant Hamiltonian approach in terms of multisymplectic frameworks may be found 
miiii- 

In mechanics one can construct solutions of the canonical eqs. of motion if one has an 
appropriate solution of the corresponding Hamilton-Jacobi equation, the solutions of 
which describe wave fronts which are transversal to a "field" of extremals and which 
contain the same dynamical information as the extremals themselves. 

For field theories this is no longer true fflj: Solutions of the Hamilton-Jacobi equation 
(HJE) associated with one of the canonical frameworks mentioned above provide "ve- 
locity" fields which in general do not obey the necessary integrability conditions (IC). 
The latter have to be postulated separately and give rise to an additional set of partial 
differential equations. However, if these equations and the associated Hamilton-Jacobi 
equation are satisfied, then combined they contain the same dynamical content as the 
Euler-Lagrange equations [p]]. Thus, it is possible to construct solutions of the Euler- 
Lagrange field equations (ELE) by solving a Hamilton-Jacobi equation and a set of 
integrability conditions. 

The aims of the present paper are the following ones: 



1. All the different canonical formalisms just mentioned coincide if the field has just 
one real component |L), especially for fields in 2 spacetime or Euclidean dimensions. 
In this case one can use the more simple DeDonder & Weyl HJE in order to 
solve it simultaneously with one integrability condition. Here we do not construct 
completely new solutions of the ELE but find one-parameter extensions of a given 
solution ip . The method is applied to the following one-component models in 2 
dimensions: massive and massless Klein-Gordon, Sinh-, Sine-Gordon, Liouville 
and 4 theory. The solutions of the HJE plus IC are constructed by expanding 
the solutions of these equations in powers of the field variables. This leads to 
a hierarchy of nonlinear PDEs that can be transformed into linear PDEs with 
nonconstant coefficients. By applying integrable Backlund transformations these 
PDEs can be reduced further to free (massive or massless) Klein-Gordon equations. 
Remarkably, it is only necessary to solve just two linear PDEs in order to obtain 
the general solution for every order of the hierarchy! 

2. The solutions of the DeDonder & Weyl HJE do not have the appropriate transver- 
sality properties required to construct the wave fronts associated with a given one- 
parameter set of extremals. Those have to be given in terms of the solutions of the 
more complicated HJE of Caratheodory. As we are dealing with one-component 
fields here only, the two canonical frameworks are equivalent and one can construct 
the 1-dimensional Caratheodory wave fronts from the solutions of the DeDonder 
& Weyl HJE. This procedure will be outlined and applied to given solutions (e.g. 
solitons) of models mentioned above. 



2 



The paper is organized as follows: 



In chapter 2 we very briefly summarize Lepage's reformulation of canonical mechanics. 

In chapter 3 we apply this framework to field theories with one component in two di- 
mensions and present the HJEs of DeDonder & Weyl and of Caratheodory, respectively. 

In chapter 4 we study the hierarchy of equations derived within the Hamilton-Jacobi 
framework of DeDonder & Weyl and the associated integrability condition for the velocity 
fields by expanding the wave fronts in powers of the field variable (p in the neighbourhood 
of v?o- 

In chapter 5 we determine the Backlund transformations which reduce this hierarchy of 
linear PDEs with nonconstant coefficients to free field equations. 

In chapter 6 we apply the general results to several well known models and point out 
some relations to stability problems of solitons. 

In chapter 7 we discuss the relations between the hierarchies of PDEs derived from 
expanding the HJE and the IC and a corresponding expansion of the Euler-Lagrange 
equations. Solutions of the equations of motions are determined from those of the HJE 
and the IC perturbatively. The kink solution of the Sine-Gordon equations is treated in 
considerable detail here. 

In chapter 8 we calculate the wave fronts by means of Caratheodory's formulation, es- 
pecially for the kink solution just mentioned. 



2 Lepage's Canonical Formulation of Mechanics 

We very briefly recall Lepage's main idea of introducing the canonical formalism in 
mechanics for one configuration variable q. The general case is discussed in [ED]. 

The starting point of the canonical theory for Lagrangian canonical systems is an action 
functional. In mechanics it is given by a Lagrangian 1-form u integrated along a path 
C:={t,q = q(t)} in the extended configuration space Aii+i := {(t,q)}: 



As to the variational principle it is preferable to consider the generalized velocity v as an 
independent variable, which coincides with q on the extremals C = C := {(t,q = q (t))}. 
Normally this extension in the number of variables is performed by using Lagrangian 
multipliers. Lepage's reformulation of the variational principle is similar in spirit Q. 
The initial canonical Lagrangian form u = L(t,q,q) dt is extended by the product of 
a Lagrangian multiplier h(t, q, v) and the Pfaffian form g = dq — v dt vanishing on the 
tangent vectors of the extremals, which ensures the identification of v(t) with q(t) on 
the solutions of the equation of motion. Then the action integral A[C] over the path 
C:={(t,q(t))} can be modified: 




(2.1) 




(2.2) 



without changing the Euler-Lagrange equations and their solutions C = C - 
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The form g generates an ideal I[g] in the algebra A of forms on the extended configuration 
space A4 1+1 := {(t,q)}: if a G A and (3 G I[g], then a A (3 is also an element of the ideal 



The Lagrangian multiplier h(t, q, v) can be fixed by varying the action integral (2.2) 
with respect to q, v independently. This leads to the standard definition of the canonical 
momentum: h = d v L=:p. 

We obtain the same results by requiring dQ to be an element of the ideal I[g\: 

dtt = (d v L - h) dv A dt + (dh - d q L dt) A g = (d v L - h) dv A dt + (modI[g\) . (2.3) 

= 0! 

Hence dfl is a closed two form on families of extremals covering the extended con- 
figuration space A^i+i = {{t,q)} - or correspondingly - a (Lagrangian) submanifold 
Q := {t, q,p — ip(t, q)\(t, q) G M} of the extended phase space P2+1 : = {(t,Q,p)}- 

Following Poincare's lemma Q is locally (at least) exact Q = dS(t,q). 

The Legendre transformation L — > H , v — > p can be implemented as a change of basis 
in the cotangent bundle T*(A4 1+ i), g —> dq, dt — > dt: 

Q = L dt + p g = L dt + p (dq - v dt) = -(pv - L)dt + pdq = -H dt + pdq. (2.4) 

H denotes the usual Hamiltonian H = pv — L = H(p,q,t). 

The existence of a potential S(t, q) for the basic differential form Q = dS yields the 
familiar Hamilton-Jacobi equation for S(t, q) and the corresponding condition for the 
momentum: 

tt = -H (t, q,p = ip(t, q))dt + ifj(t, q) dq = dS(t, q) = d t S(t, q) dt + d q S(t, q) dq . (2.5) 
Comparing the coefficients of dt, dq yields: 

d t S(t, q) + H (t, q,p = i/>(t, q)) = 0, P = ip(t, q) = d q S(t, q) . (2.6) 

The extremals can be determined, if a complete integral of the Hamilton-Jacobi equation 
is found. Dealing with one independent variable q this integral depends on one constant 
cq- The solution of the equation of motion can be calculated by solving the algebraic 
relation d CQ S(t, q, cq) =c\ with a second constant c\. 

In mechanics there exists a special relation between the wavefront S(t, q) = a = const, o G 
TZ and the extremals q = qo(t) called "transversality": the union of their respective 
tangent spaces span{e t = d t + q(t)d q } and span{w=pd t + Hd q }, span the whole tangent 
space 7p(A^i+i) at any point P G M1+1, if the Lagrangian L does not vanish, because 

of 

1 p 
v H 



det (e t , w) =det 



{H-qp) = -L. 



This concept of the Hamilton-Jacobi framework developed for mechanics can easily be 
generalized to field theories. We confine our discussion to those theories depending on 
one real scalar field tp = tp(z, z) in a 1 + 1 dimensions. Here z, z play the role of lightcone 
variables. For details see ffl ITOfl [BH. 



3 The Hamilton— Jacobi theories of DeDonder & 
Weyl and of Caratheodory 



As in mechanics a canonical theory for fields is based on an action functional A[T] 
defined on a two dimensional surface E:={(z, z, tp(z, z))} in the extended configuration 
space M2+1 -={(z, z, <p)}: 

A[T] = J to = J C(z,z,ip(z,z),v = d z <p,v = d z ip) dz Adz . (3.1) 

Only on the extremals <p = ip (z,z) the generalized velocities v,v coincide with the 
derivatives of the fields: v = d z (po,v = dzipo. 

The canonical 2-form uj = £ dz A dz is extended by means of two Lagrangian parameters 
h(z, z, <p), h(z, z, ip) and a 1— form g = dip — v dz — v dz that vanishes on the 2-dimensional 
extremals ip = (p Q (z, z): 

Q = Cdz /\dz + hdz A g + hg Adz . (3.2) 
The Lagrangian multipliers h, h are determined by requiring dQ G I[g]: 

dQ = {d v C ~ h)dv Adz Adz + (d d C -h)dv Adz Adz + {modI[g]) = (modI[g\) , (3.3) 
yielding 

h = d v C =: p , h = d, d C =: p . (3.4) 

The Legendre transformation C — >• 7i, {v,v} — > {p, p} can be implemented as a change 
of the basis in the cotangent bundle T*(A4 2 +i), Q dip, dz dz, dz — > dz: 

Q = — 7i dz A dz + pdz A dip + pdf A dz with TC := pv + pv — C. (3.5) 

Because dfl = (modJ[^]) it is locally exact Q = dS, S G T*(A^2+i) on families of 
extremals. However, contrary to mechanics Q as an exact 2-form can be represented in 
different ways by means of a Pfaffian form S. In the case of DeDonder & Weyl [Kj |TT[ : 

<Sdw — S(z, z, ip) dz — S(z, z, ip) dz (3.6) 

and in the case of Caratheodory |TI|: 

S c = S'(z,z,ip)dS*{z,z,<p). (3.7) 

Comparing the exterior derivatives of these expressions with equation ( |3.5| ) 

Q = -Hdz Adz + pdz Adip + pdip Adz = dS A dz — dS A dz , (3.8) 

we obtain the Hamilton-Jacobi equations and the transversality conditions for a one 
component field theory in DeDonder and Weyl's formulation: 

d z S + d- z S = -H , P = d^S, p = 3^S (3.9) 

and in Caratheodory's case: 

tt = —H dz A dz + p dz A dip + p dip A dz = dS z A dS 2 , (3.10) 
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we get: 

n — f)- ^ z f) <? 2 — f)- ^ z f) ^ z 

o z b a- z b a z b o- z b - n, p = dzS z d(pS z _ dzS z d<pS z • ^.nj 

The two theories here are equivalent because a n-form in a space of n + 1 variables has 
always rank n ||. Due to this algebraic equivalence of covariant canonical formula- 
tions for one component field theories we may choose DeDonder and Weyl's description 
to embed the extremals of interest in a system of solutions of the Hamilton- Jacobi 
equation. We use Caratheodory's theory for the explicit calculation of the wave fronts 
S z = const., S z = const.: 

In two dimensional field theories involving one field variable the basic two form Q has 
always the rank two, i.e. it can be constructed from two independent one forms by lin- 
ear combination of exterior products. Because Q is closed its rank is equal to its class, 
that gives the codimension of the integral submanifold determined by Q. To calculate 
this integral manifold — the wavefronts in our case — we can use a corollary of Frobe- 
nius' integrability theorem |IJ: there exist two functions S z (z,z, <p), S z (z,z,{p) such a 
manner that the exterior product of their differentials equals Q. The corresponding 
one dimensional wave fronts are just the submanifolds determined by the conditions 
S z (z, z, ip) = const., S z (z, z, ip) = const.. 

Thus the wave fronts are equipotential surfaces of the solutions of the Hamilton- Jacobi 
equation formulated in Caratheodory's framework. For simplicity we first solve the HJE 
of DeDonder & Weyl and the associated IC and afterwards we return to Caratheodory's 
equation in order to obtain an explicit expression for the wave fronts. 

In mechanics for one variable q it is possible to construct wave fronts for 1-parametric 
families of extremals that cover a certain region of the extended configuration space Aii+x 
and vice versa. Provided a solution S(t, q) of the Hamilton-Jacobi equation (HJE) is 
given, the corresponding velocity field, the so-called " slope function" : 

$(t,q) = d p H(t,q,p = d q S(t,q)) (3.12) 

determines the corresponding 1-parametric extremals by means of the ordinary first 
order differential equation: q(t) = &(t,q(t)). 

In general this is not true for field theories; the ability to embed extremals lpq(z,z) 
in a given wave front can be maintained only if the slope functions (velocity fields) 
v = <$>(ip, z,z), v = <&(ip,z,z) obtained from the inverse Legendre transformation 

d z ip(z,z) = v (p = d v S,p = dpS^, z, ipj = z,z) , 

d z (p(z,z) = v (p = d v S,p = d v S,z,z,ipJ = $(ip,z,z) (3.13) 

satisfy the integrability condition d z d 2 (p(z, z) =d z d z (p(z, z), which results in the require- 
ment 



-^$(z, z, <p(z, z)) := d z § + $ • <9 V $ = z, <p(z, z)) := d z <£> + $ • <9^$ 

az dz 



(3.14) 



on $ and $. 

Provided the Hamilton-Jacobi equation and the integrability condition are fulfilled, the 
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Euler-Lagrange equation is satisfied automatically. This can be seen as follows: 

Differentiating the Hamilton-Jacobi equation ( |3.9| ) with respect to the field variable ip 
one obtains: 



d z p + d z p 



-d v Tt 



Because d^K- 



-d^C, d p H- 



$ and d p H 



- d p Hd v p - dpHd^p . 
$ (see [0) we get: 



d d _ ddC d DC 
—p + — p = — — + 

dz dz 



(3.15) 



(3.16) 



dz dv dz dv 

The momenta p(z, z, (p) and p(z, z, (p) are defined on the extended configuration space 
M.2+1 an d are considered to be associated with a family of extremals if = (p(z,z,u) 
parametrized by u. The Euler-Lagrange equation has to be fulfilled for every single 
extremal tp={p(z, z, u = const.). Hence if we insert p and p into this equation of motion, 
we have to be aware of their implicit dependence of z, z via the field variable (p. Taking 
the defining equations ( |3.13|) for the slope functions into account, the total derivatives 
in the PDEs ( p,14| ) and ( |3.16| ) are defined as: 



d 

dz 



d_ 

dz 



=const 



dtp 



=const 



d 

di 



d_ 

dz 



=== d 

=const at P 



(3.17) 



=const 



Therefore the total derivatives are nothing but derivatives with respect to the indepen- 
dent variables z, z regarding u to be a constant. 

Due to the extend in which the integrability condition is taken into account there exist 
two different methods of using a Hamilton-Jacobi theory: the weak and the strong em- 
bedding of extremals in families of wave fronts. 

1. ) Weak embedding: 

This method is used to embed a given single extremal <p>o{z, z) in families of wave fronts. 
In order to obtain a weak embedding it is sufficient to take only the Hamilton-Jacobi 
equation and the transversality conditions (|3.13|) on Q (z,z) into account. In this case 
usually one chooses a linear ansatz in the field variable for one of the functions S or S. 
However, this approach in general will not provide new extremals, because the IC are 
fulfilled on the given extremal only. For details to this subject see [|T9|]. 

2. ) Strong embedding: 

Here one requires the IC not only to hold on the given extremal but in a whole neigh- 
bourhood of it. If this is the case then one is able to generate new extremals from a 
given one by integrating the slope functions ( |3.13| ). 

In the following we study the strong embedding of a single given extremal into a family 
of wave fronts. 

Like in mechanics there exists a transversality relation between extremals (p(z, z) and 
wave fronts S z (z,z,ip) = const., S z (z,z,(p) = const., which holds, if in every point 
P G M. i+2 the basis of tangent space TM. 1+2 is given by a union of the basis of the 
2-dimensional tangent space of the extremals e 2 = d z + v d v , ez = d 2 + vd^ and a basis 
vector w =pd z + pd z + Tid^ of the 1-dimensional tangent space of the wave fronts, i.e. 
iff the Lagrangian density L does not vanish: 



det 



= det 




(TC — pv — pv) = —C 7^ . 



(3.18) 



Notice that in theories with more than one real field (d > 2) both £ and 7i have to be 
nonvanishing quantities to guarantee the transversality relation 



4 Hamilton— Jacobi theory for one real field in De- 
Donder and Weyl's formulation 

We here restrict ourselves to Lagrangian densities of the following type: £ = d z ipd z ip — 
V(f). The potential V(tp) is an analytic function. Here we have the canonical momenta 
p = v, p = v , the Hamiltonian density H = pp + V and the slope functions $ = d v S, 
§> = dpS. We have the (DeDonder & Weyl) Hamilton-Jacobi equation 

d z S + d- z S = d v S d v S + V{ip) (4.1) 

and the related integrability condition 

a a d v s - d- z d 9 s = o v sdls - d v s%s . (4.2) 

Knowing solutions S and S of the equations ( |4.1|) and ( |4.2j ) a family of embedded ex- 
tremals ip = (p{z, z) is determined by a system of first order PDEs: 

d z ip(z,z) = $ = d lp S{z,z 1 (p = (p) , (4.3) 
dg<p(z, z) = $ = d v S(z, z,<p = <p). (4.4) 

A solution is obtained by expanding S(z, z, ip) and S(z, z, <p) in powers of the difference 
y = ip — tpo between ip and a known extremal <po(z, z): 

OO J OO 1 

S(z,z,<p) = ^A n (z,z)y n , S(z,z,ip) = Y / ^M^z)y n . (4.5) 

This method of expanding about a given solution of the equations of motion is commonly 
employed e.g. with stability investigations or determining (quantum) fluctuations around 
c-number fields in selfinteracting theories |13|] |T5 . 



Naturally <po has to satisfy the transversality relations (P), that fix only the functions 



A\ = d z (po and A\ = d z <po, without influencing the remaining coefficients. Inserting the 



expressions (P) into the HJE ( |4.1| ) and expanding the potential V in powers of y we 
see that the equation is automatically fulfilled up to the order y 1 , whereas the IC ( | . 2| ) 
is fulfilled on the extremals up to the order y°. 

A and A Q are only affected by the HJE of zeroth order in y: 

d z A + d z A = C\ v=(po = d z (p d z ip - V(<p Q ). (4.6) 

The general solution of the homogeneous equation d z AQ° m + d z A^ om = is given — at 
least locally — by Aq 0171 = d 2 Xo( z , z) an d A^ " 1 = —d z Xo(z, z). One special solution of 
the inhomogeneous PDE ( f4.6|) can be obtained by integration: A™ h = J C\ Q dz/2 and 
A l Q h = J C\ o dz/2. This yields the general solution (locally): 

A (z, z) = ^ y [d z <p d z ^ - V(ip )} dz + d z xo(z, z) , (4.7) 

A (z,z) = -J[d z <p d z (p -V(<po)]dz-d z Xo(z,z). (4.8) 
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In order to determine the coefficients A 2 and A 2 , we insert A\ = dz<fio and A\ = d z ipo into 
the IC of order y: 



d z A 2 - A 3 d z (f + AiA 3 + A 2 A 2 = d- z A 2 - A 3 d z <p + A 1 A 3 + A 2 A 2 



(4.9) 



Since y is a function of z, z, its derivative with respect to z or z yields d z y = —d z ip 
and d z y = —d z ifo respectively. Thus we infer from equation fl4.9|) that d z A 2 = d z A 2 
which permits to express these two functions at least locally by one generating potential 
function: A 2 = d z ( 2 (z, z) and A 2 = d z ( 2 (z, z). ( 2 has to fulfil the PDE 



d z d z ( 2 (z, z) + d z ( 2 (z, z)d z ( 2 (z, + ^ [d?V ((p=(p (z, z)) 

The ansatz 

C 2 = \n9(z,z) A 2 = d z \n9, A 2 = d z \n6 

linearizes the HJE of second order in y: 







d- z d z 9(z,z) + - &,V(<p = <p (z,z)) 6(z,z)=0. 



(4.10) 
(4.11) 

(4.12) 



In order to obtain the expressions for A n and A n one has to substitute the power series 
( fl.Sp into the IC and the HJE and to compare the coefficients of the powers y n ~ l and 
y n respectively. Starting with the IC in (n — l)-th order: 

n-l 



d z A n — A n+ i 

p=0 

n-l 

d z A n - A n+1 d z (p + 

p=0 



n — 1 

V 

n — 1 
P 



Ap+iA n - p+ i 



(4.13) 



we draw our attention to the coefficients A n+ i and A n+ i of highest order. Fortunately 
they disappear from this equation as well as from the HJE due to the relations A\ = d z (po 
and A\ = d z (fo. Provided the coefficients A , . . . , A n _i and A , . . . , A n _i have already 
been determined recursively, one gets an equation for the functions A n (z, z) and A n (z, z): 



n-2 



d z (e n - 2 A n )-d- z (9 n - 2 A n ) = 9 n - 2 Y 



n — 1 

. V . 



n — 1 
p-1 



Ap+lAn. 



■p+1 



p=2 L 

=: lnh(A 0} ... } A n _ 1 ;A 0} ...,A n A . 



(4.14) 

The Inhomogeneity Inh (^A 0} . . . , A n ^i, A 0} . . . , A n _-^j vanishes for n < 4. Similarly to 
( |f .8|) this hierarchy of equations is solved by splitting the solution in two parts: a general 
solution of the homogeneous part d z (6 n ~ 2 A^ om ^j—d z (6 n ~ 2 A^ om ^j =0, that is obtained - 
at least locally — by 



.4 



horn 



l ' n d z Cn{z, z) and A 



horn 



6 2 - n d z ( n (z,z) 



(4.15) 



with an arbitrary smooth function ( n (z,z). The second part is a special solution of the 
inhomogeneous equation ( (4.14| ): 

(4.16) 
(4.17) 



A 



inJi 



6 2 - n d z Xn(z,z) and AT = -e 2 - n d z Xn(z,z) with 



Xn{z,z) = - 



Inh [A Q , . . . , A^i, A , 
9 



A n _i) dz' dz' . 



The general solution 



A n = d 2 - n d- z [( n (z,z) + Xn(z,z)] and A n = B 2 ~ n d z [C n (z, z) - x„(z, z)] (4.18 



has to satisfy the HJE ( |4.1[ ) to n-th order in y, in which, remarkably, the coefficients of 
highest order A n+ i and A n+ i vanish as well as in the IC, due to the relations Ai = d 2 fo 
and Ai = d z ifQ. Separating the functions of the remaining highest order A n , A n leads to 
the equation: 



n-2 



d z (6 n A n ) + d- z (e n A n ) = -e n J2 ( " ) A p+1 A n _ p+1 - 9 n -j-^V {<p) 



p=2 



p, 



d{p % 



(4.19) 



tp=ipo 



It is convenient to set ( n (z, z) = Q~ x Xn before inserting into expression (|4.19| ). Using the 
relation ( [4.12| ) the HJE of n-th order in y yields an equation, which is nothing but the 
inhomogeneous extension of the linear PDE ( |4.12| ): 



d z d z \n + ^ d 2 V{ipQ 



Xr, 



Inh (xa,...,Xn-i) 



(4.20) 



with the inhomogeneity: 



n-2 



Inh ( Xo , • • • , Xn-i) = d- z 6d zXn - d z 6d- z xn - ^9 n ~ l 

A p=2 



n 



Ap+iA n _ p 



+i 



(4.21) 

In general the coefficients A n , A n are determined by the n-th order of the HJE and the 
n—1 order of the IC for n > 3: 



A n = e 2 ~ n d- Z 



Xn{z,z) 
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Xn{Z,Z) 



A n = e 2 - n d, 



Xn{z, Z) 

9 



+ Xn{Z,z) 



(4.22) 



Notably the infinite hierarchy of functions Xn{z, z) has to fulfil only one PDE of second 
order: the equation ( |4.20| ). Its integral can be obtained by determining the general 
solution of the homogeneous PDE ( J4.12 ), which is the same for all orders n > 2 and 
one special solution of the inhomogeneous equation ( 4.2C ). It is given by using a Green 
function that can be chosen to be the same for all orders n without any loss of generality. 
The solutions of (4.12) can be obtained by employing Backlund transformations. 



5 Backlund transformations 



Backlund transformations (BTs) are maps between the tangent bundles of integral sub- 
manifolds associated with PDEs. If we are able to find a BT from the PDE which we 
wish to solve and to another one the general integral of which is known, we can obtain 
the general solution of the first one by integrating the BT. This treatment of a single 
PDE can be generalized to systems of partial differential equations |§. 

By applying BTs we want to reduce the linear PDEs of second order with nonconstant 
coefficients of type ( 4.12 ) to linear PDEs with constant coefficients. The inhomogeneous 
extension of equation (|4.12| ) could be solved by BTs, too, but for sake of simplicity we 
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construct a special solution of the inhomogeneous equation (|4.12|) by using a Green func- 
tion and Fourier transformation. Like other authors, e.g. f?J, the general form of a BT 
we are starting from is given by two functions F\, F 2 : 



d z 9(z,z) = Fi[z,z,9(z,z),9(z,z),d z 9(z,z) 
d z 8(z,z) = F 2 z } z,8(z } z),8(z } z),d z 8(z,z) 



(5-1) 



9 has to fulfil the relation (|4.12|) , whereas 9 denotes the transformed function which is 
supposed to obey a linear PDE with a constant coefficient m 2 : d z d z 8 — m 2 8 = 0. Of 
course m 2 can be equal to zero which would yield the wave equation. Obviously F\ and 
F 2 have to fulfil the integrability condition d z d z 8 = d z d z 8 =>- dF^jdz = d_F 2 /dz. This 
integrability condition does not lead to a restriction on the solutions of the PDE ( [4.12 ), 
if such a BT is found. Thus Fi and F 2 have to fulfil the two equations 

diq/df = dF 2 /dz , dF^dz = m 2 9 . (5.2) 

We would like to point out that the more general ansatz d z d z 9 = K(z, z)9 leads to the 
same BTs Q5.9Q below at least for the models discussed here |T8|| . 



Differentiating the eqs. (5T2") with respect to d z 9 and d z 9 leads to 

d 2 F = , =}► F 1 = f 1 (6, 9, z, z)d z 9 + mi(6, 9, z, z) 



d 2 F 2 



0, 



f 2 (9, 9,z,z)d- z 9 + m 2 (9,9,z,z) 



(5.3) 
(5.4) 



Substituting these expressions into eqs. ( |5.2|) , comparing the coefficients of d z 9, d z 9 
and differentiating with respect to 9, 8, we conclude by lengthy but straightforward 
calculations that 



fi = ~ J 2 = c , c eTZ, c ^ , 

= dgmi-fidgrrn, mi = gi(fi9 + 9, z, z) , (5.5) 

= d e m 2 + f\d~ e m 2 , m 2 = g 2 (fi9 - 9, z, z) . (5.6) 

Inserting these results into relations ( |5.2j ) and differentiating with respect to fi9+9 = i]i 
and f\9 — 9 = r\ 2 we obtain: d 2 l m 1 = 0, d 2 2 m 2 = 0. Substituting f x 8^8, j\d z 8^d z 8 and 

f\d- z 9^d z 9 we get: 



d z 9 
d-J 



d z 9 + ai{z,z) 9 + 9 +p 1 (z,z) 



—d z 9 + a 2 (z, z) 



9-8 



+ P2(z,z) 



(5.7) 
(5.8) 



Considering the two functions /3i, f3 2 is only necessary for transformations between in- 
homogeneous PDEs. Thus, we here can choose (3i = = (3 2 . Inserting our results ( |5. 7| ) 
and ( |5.8|) into equations ( |5.2|) and comparing the coefficients of 9, 9 we finally obtain the 
linear BT 

' ' ' ' (5-9) 



dj = d z 8 + {& + §} d z ip , d z 9 = - d- z 9 + {§ - 8) d z ip 



with the BT generating function ip = ip(z, z) which is a special solution of 

1 



d z d- z ^j - (d z ^) («9#) - -djy&o) = , dM + (d z ?P) (d z iP) - m l 
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(5.10) 



Substituting ijj = — In 



in these equations we conclude 

r 



dM + wv> = o , o z d- z ^ 



(5-H) 



We therefore have to solve the following problem: 

2mm We /iai>e to find a solution of the equation (|4 . 1 2|) the inverse of which has to fulfil 
a Klein-Gordon or a wave equation, then we can integrate the linear BT and obtain the 
general solution of ( [4.12 ). 



Another method of solving the eqs. (|5.10|) is based on their linear combinations: 

d M drf> = \ (2m 2 + flj%)) , (d^)(d-^) = 1 (2m 2 - 3}%)) . (5.12) 

Integrating the first equation and inserting this solution into the second expression im- 
poses a restriction on d^V(ipo)- We make use of eqs. ( |5.12j ) when we study the 4 -model. 



6 Applications 



After solving the homogeneous equations ( f4.12|) , we calculate a special solution of their 



inhomogeneous extension ( 4.20|) by determining a Green function — without need of 



specifying the inhomogeneity Inh. To obtain a solution of equation ( |4.20|) specific for the 
models under consideration we have to fold the Green function with the inhomogeneity 
in every order y n . 

The Hamilton-Jacobi theory for the non-selfinteracting scalar field theories Co = d z (pd z f 
and Ci = d z (pd z ip — l/2m 2 (p 2 with the light cone variables z = (x + t)/2 and z — (x — 1)/2 
leads to the wave or the Klein-Gordon equation (|4.12|) without need for a Backlund 
transformation or specifying an extremal <po(z, z): 

£o'- d z d z Xn = Inh (xo, • • • , Xn-i) , C\: d z d z \n + m 2 Xn = Inh (xo, • • • , Xn-i) ■ (6-1) 

The general solutions of these relations are known. Therefore we draw our attention to 
the more interesting case of selfinteracting theories: 



6.1 The homogeneous equations 
6.1.1 Liouville model 

Applying our formalism to Liouville's theory £ = d z (pd z (p + 2 exp((p) , using an arbitrary 
solution of the equation of motion for which the general expression is known 0: 

^{ 2 (M^}, ( , 2) 

with arbitrary smooth functions s(z) and s(z), the relation ( |4.12| ) yields: 

d z d 2 9 L -2 ^ zS) } d jf ) e L = 0, d s d- s 6 L (s,s)-2—^d L (s,s) = (6.3) 
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by employing a transformation of variables z — > s(z),z — > s(z). Obviously one special 
solution of this equation is 6$ = l/(s+s). Its inverse fulfils the wave equation d s d s 9 L = 0. 
Thus we know that there exists at least one BT which connects the integral submanifolds 
of (|6.3|) and of the wave equation. Returning to the equations (|5.11|) we conclude that 



ip L = \n(s+s). So we can determine the BT by integrating the relations ( |5.9| ) and obtain 
the general solution of the equation (|6.3|): 



6 L (s, 7s) = d s C(s) - d s C(s) + __-{C(a) - C(s)} (6.4) 

[s + s) z 

with two arbitrary smooth functions C(s) and C(s). 



6.1.2 Sine & Sinh— Gordon model 



In two dimensional space-time two types of solitons exist: the "bell" with the same 
asymptotic value at x = — oo and x = oo and a "kink" soliton with different asymptotic 
values. Moreover there exist a topological conserved quantum number associated with 
the asymptotic behaviour of these solitons. The corresponding conserved current is given 
by: J M = e^dyip with the antisymmetric tensor & v = — e" M , fi, v — 0, 1. Thus the charge 

/oo 
J da; = (/9| a ; =oo — <^|a;=-oo) which vanishes obviously 
-oo 

in the case of the bell solitons. For kinks it is a non-trivial quantum number. 

The Sine-Gordon theory possesses a infinite hierarchy of multikink solutions, which can 
be constructed by Auto-Backlund transformations. In this model the conserved quantity 



is associated with a particle number. For details as to solitons see e.g. Ul3| , flH |, and 



0- 



We want to embed for the Sine-Gordon model Li = d z (pdz(p + 2[l — cos((p)} the (anti- 
) kink solution (p = ±4arctan[exp(z + ^)] and in the case of the Sinh-Gordon model 
£3 = d z (pdz(p — 2[1 — cosh(</3)] the bell solution (p = ±4 arctanh[exp(;2 + z)], which is only 
defined for z + z <0. We then obtain for equation (|4.12|) : 

Sine-Gordon: d z d- z SG - {2 tanh 2 (z + z)- 1}6 SG = , (6.5) 

Sinh-Gordon: d z d z 6 Sh - {2 coth 2 (z + z) - l}6 Sh = , (6.6) 

respectively. Following the discussion of Liouville's theory we are able to solve these two 
equations by one BT. For this it is sufficient to realize that the inverse of the two solutions 
9q G = 1/ cosh(z + z) and 6q H = 1/ sinh(z + z) solve the Klein-Gordon equation d z d z 6 = 6. 
The functions 6q G = if) SG and Oq 11 = tp Sh can be calculated by using the relations ip = 
exp(— ip) and ( |A.21| ), ( A.2 ) shown in the Appendix. Thus the two generating functions 



are ip SG = \ia(cosh.(z+z)) and ip Sh = \n(smh(z+z)) which determine the BTs (|5.9| ) between 
the Klein-Gordon eq. and the eqs. (|6.5|), (|6.6|). Their general solutions can be calculated 
by integration of the linear BTs ( |5.9| ): 

/oo roo I 2 2i(7 1 

/ exp[-i(qz + qz)}5(qq + 1)T SG (q,q) \l — — - — -f- tanh(z+z) dqdq 
-00J— 00 q + 1 q + 1 



+c cosh (z + z) , (6.7) 
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e sh = f°° r exp [_i( kz + kz)]S(kk + l)T Sh (k, k) K fc + fc ) _ cot h(z + z)\ dk dk 

+cisinh _1 (z + z) , (6.8) 

with arbitrary constants Cq,Ci G 1Z and two arbitrary functions T SG (q,q),T Sh (k,k) 
which have to be chosen in such a way that the integrals exist. A property of this BTs 
is that the solution which was used for the transformation, is multiplied with a constant 
and added to the modified solution of the Klein-Gordon or the wave equation. 

This static kink y?o = ±4 arctan[exp(x)] can be transformed by a Lorentz boost into 
time dependent solutions of the relativistic invariant Euler-Lagrange equation ipo = 
±4 arctan[exp(7(x— vt) + 5)}, 7 2 = 1/(1— v 2 ) parametrized by the velocity v and the phase 
shift 5. We are able to include the embedding of these solutions in our discussion making 
use of the transformation of variables 

z -> w = zj(l - v) +5/2, z -> w = zj(l + v) + 5/2 . (6.9) 

The same holds for the static solitons of the relativistic covariant Sinh-, and 4 -models. 



6.1.3 4 -model 

Contrary to the three previous models the following ones can only be solved by at least 
two BTs: the 4 -theories 

I) C = d s (pdsfp - 4(^ 2 + 2tp 4 and II) C = d z tpd s (p + 2tp 2 - 2^ . (6.10) 

Both theories have soliton solutions. We choose the (anti-) kink ip I = ±tanh.(z + z) and 
in case II the bell solution cp^ 1 = ± cosh" 1 (2 + 2). Then we get the two relations ( |4.12[ ) 
for 9 1 and 9 H : 

I) : d t d g 9 J - {6tanh 2 (z + ^) - 2}6 I = 0, (6.11) 
II) : d z d- z 6 H - {Qtanh 2 (z + z) -5}6 H = 0. (6.12) 

Except for the special values of same of the constants they are the same PDEs as that 
of the Sine-Gordon model, but they cannot be solved by one BT only (see Appendix). 
Therefore we employ two BTs for each model: the first BTs leads to two eqs. in which the 
coefficient in front of tanh 2 (z+z) is reduced to 2, the same as the one in the Sine-Gordon 
theory. This allows us to obtain two Klein-Gordon eqs. which differ by the choice of m 2 
after a second BT for each model. 



Inserting the function d 2 V((fo) = v (z+z) in the results (|A.29|) and ( A.2|) of the Appendix 



we are able to calculate the functions ip 1 , ip 11 and the two solutions of the equations 
( p-llQ , (p. 12 ) which are necessary for the transformations: 



I): 0o = cosh" 2 (z + z) , II) el 1 = exp[d {z - z)\ cosh" 2 (z + z) , d 2 = 3 . (6.13) 

Thus, with ip 1 = 2\n{cosh(z + z)} and ip 11 = d (z — z) + 2 ln{cosh(z + z)} we obtain after 
one BT and denoting the function 6 of eq. Q5.9Q by 9: 



I) : - {2 tanh 2 (z + z) + 2} 9 1 = 0, (6.14) 

II) : d z d- z 9 H - {2tanh 2 (z + ^) - l}9 n = 0. (6.15) 
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The second relation is the same as in the Sine-Gordon theory. Thus we only have to treat 
the first case here. This PDE has the special solution 9q = 1/ cosh(z+z) exp[idi(z— z)], d\ = 
3 the inverse of which fulfils a Klein-Gordon equation, namely: d z d s (l/9l) = 4/0g. 
The solution 9q was found according to the method discussed in the Appendix. The 



generating function for this BT is = ha{cosh.(z + z)}+id 1 (z — z) 
resulting solutions of the eqs. ( |6.11| ) and (|6.12|) are: 



\tl(9q). Hence the 



9 1 



/oo roo , 
/ T% q)6(qq + 4) exp[-i(qz + qz)] { q 4 + 16 - 4q 2 
-oo J — oo 

- (6iq 3 - 24iq) tanh(^ + z) — 12g 2 tanh 2 (z + z)} dqdq 
+co cosh~ 2 (z + z) + c\ cosh - (z + z) tanh(z + z) exp[id\{z — z)\ , (6.16) 
/ T H (k, k)S(kk + 1) exp[-i(kz + kz)] \-{k 2 + l) 2 

-oo J —oo ^ 

-6iJfc(l - A; 2 )tanh(^ + z) + 12A; 2 tanh 2 (z + z)\ dkdk 

+c 2 cosh~ 2 (z + z) exp[d (z — z)\ + c 3 cosh _1 (z + z) tanh(z + z) (6.17) 

with d 2 ] = dl = 3; Co, c±, C2, 03 <E7Z and two arbitrary functions T / (g, q), T n (k, k). 
Expressions like ( |6.5|) and (|6.11|) also occur e.g. in stability investigations or in discus- 



sion small fluctuations around the known soliton solutions y?o of these theories [13[ fli~5 



Carrying out the second variation of the action functional in the case of the Sine-Gordon 
and the 4 -theory employing the soliton solutions given above — or equivalently insert- 
ing (p = ip +ip c , \ip c \ <C \<p \ into the Euler-Lagrange equation of the Sine-Gordon or the 
4 -theory — yields: 



— <9 2 + n(n + 1) tanh 2 (x) + Co u{x) = ui 2 u(x) 



(6.18) 



where (p c (x,t) = exp(iut)u(x) . We have n — 1, cq = 2 for the Sine-Gordon and n — 2, 
Co = —2 for the 4 -model. Stability of the soliton solutions requires that all eigenvalues 
uj 2 of this Schrodinger-like equation should be non-negative, so that small perturbations 
about (p do not grow exponentially in time. The lowest eigenvalue is u 2 = and the 
corresponding solution C is the translation mode. It must be present, because of the 
translation invariance of our models under consideration. Remarkably we have used it 
for both models to construct the Backhand transformations: (p c —\j cosh(x = ,2+^) and 



0i- 



1 / cosh 



x 



The eqs. (|6.5| ) are (|6.11| ) reduced to eqs. (|6.18| ) if one if one sets x = z + z, t = z — z and 
exp(iut)u(x) . The eqs. ( |6.18| ) are solvable by transforming them into hypergeometric 
differential equations the solutions of which can be given by finite series in powers of 
tanh(x) functions |L7}] : 



SG: ul = k 2 + l : 



u{ x) 



exp(ikx){k + itanh(x)} 



oo z = 3 

.2 _ ;„2 



u(x) = tanh(x)/ cosh(a;) 



(6.19) 
(6.20) 



J k~ k z +4: : u(x) = exp(ikx) j — 1 — k 2 — 3Htanh(x) + 3 tanh 2 (x)|(6.21) 
These solutions are contained in our more general results (|6.7|) and ( |6.16|) . 
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6.1.4 A mathematical remark 



As discussed in the Appendix our results can be generalized in order to reduce the 
hierarchy of linear PDEs: 



dgdgO = {n(n + l)r] 2 + a}9 , aeTZ, n = 0,1,2, 



(6.22) 



by n BTs to a Klein-Gordon or a wave equation successively. One BT can raise or lower 
the coefficient n to n+1 or n—1. We have to assume that the smooth function rj(l) fulfils 
the nonlinear differential equation: dir] = brj 2 + c with c G TZ, b = ±l and l = z+z: 



- > 

c 

C 



V 



c tan 



c tanh 



c\(l + l ) 



c\(l + l ) 



cell, l eC 



: T) 



i 



(6.23) 
(6.24) 



The PDEs ( |6.22|) and ( |6.24j) are solvable without need of specifying the constants a or 



c. 



Obviously the PDEs ( |6.f 1| ), ( p. 12] ), ( |6.5| ) and (|6.6|) correspond to special choices of a, 6, c 
and n: b = —l, c=l, a = — 1 and n = l give the Sinh & the Sine-Gordon models, whereas 
n = 2, b=— 1, c=l give the 4 -theory with a = — 2 (case I) and a = — 1 (case II). 



6.2 The inhomogeneous equation 

A special solution of the inhomogeneous equation ( |4.20| ) can be obtained for the Sinh 
& Sine-Gordon and both cases of the 4 -theory by employing Fourier transformations. 
We discuss these theories first and return to the Liouville model later. 



6.2.1 Sinh & Sine-Gordon, 4 equations 



We introduce the Green function G(z, z, z, z) for the inhomogeneous equations ( |4.20| ): 

d z d- z G + \ {%V((po)} G = 6(z - z)S(z - z) . (6.25) 



If a Green function is found, we can calculate the solutions \n of eq. ( |4.20| ) to all orders 
of y n : 



oo r+oo 



Z7T J — oo J — oo 



Inh {xo(z, z),..., Xn-i(z, J)} G(z, z, z, z) dz dz . 



(6.26) 



Introducing the Fourier transform G by: 
1 



G(z, z, z, z) 



4tt 2 



oo r+oo 



exp{i[q(z — z) + q(z — z)]}G(z, z, q, q) dq dq . (6.27) 



oo J — oo 



As we need only a special G we try for G the ansatz that depends on q, q and l = z + z 
only. We then obtain 

(6.28) 



1 = dfG - i(q + q)d t G - \dlVG . 
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In our cases the potential term is only a function of /: d 2 V = 2ao+2n[n+l]f 2 (l), do G 1Z, 
n = 1,2 and / denotes a tanh(Z) (see (|6.5| ), (|6.12 ) and ( |6.11|) ) or coth(Z) (see ( |6.6D ). In 



order to obtain one solution of the PDE fl6.28|) it is sufficient to assume any / that solves 
dif = l — f 2 . Since ( |6.28 ) is only a linear differential equation we are able to calculate a 



special inhomogeneous solution, if we know a homogeneous one: 

G inh (q,q,l) = I G hom (q,q,iy 2 exp(i[q+q]l) I G hom (q, q, I) exp(-i[q+q]l) dl dl . (6.29) 

Substituting G/ lom (g, g, Z) = exp[z(g+g)Z/2]G(g, g, /) leads to: 

8?G-|ao- (9 ~^ )2 +n(n + l)/ 2 |G = 0. (6.30) 

Choosing the special ansatz: G = (c + Ci/) exp(c 3 Z) and comparing the coefficients of 
powers of / we get for n = l with a = a — (q — g) 2 /4: 



: c 3 = ±va + 2 , c = — Cic 3 , c a G 72- , Ci 7^ , (6.31) 
a = -2 : c 3 = = c , c x G ^ , c x + . (6.32) 

This choice of parameters provides the Green function for the Sinh Sz Sine-Gordon 
theories (|6.6j ) and ( |6.5|) . The case n = 2 is of interest within the 4 -models Q6.ll ) and 



( p,12| ). Here we have to add the term 02/ exp(c 3 /i) to G. This yields finally: 



: 3c = c 2 (5 + a) , c 2 = — — , c 3 = ±y/a + 6 , ci G 72-, ci 7^ , (6.33) 

c 3 

a = — 6 : c\ = c 3 = , 3c = — 2c 2 , c 2 e%, c 2 7^ . (6.34) 

The sign of c 3 has to be chosen in such a way that the Fourier integral for the Green 
function G(z, z, q, q) exists. 

Inserting these results into the expression for G m /, and inverting the Fourier transfor- 
mation we are able to obtain — in principle — the solution of equation ( |4.20| ) for every 
order of y n by integration. 



6.2.2 Liouville Theorie 



For the Liouville model we start with equation (|6.25|) , too and obtain instead of eqs. 
with relation fl6.30| ), where l — s(z)+s(z): 

{q-qf , 2 



dfG - /' + - G = . (6.35) 



4 / 



The solution is 



G = c |yO sin(aZ) — a 2 cos(a/)| + c\ jya cos(aZ) + a 2 sin(a/)| 



(q-q) 



2 

(6.36) 

Obviously we can not integrate G^ explicitly, if use of this solution. However, the 
Green function G can be obtained for the Liouville model. We return to equation ( |6.25|) , 
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choose the ansatz G(s(z), s, s(z), s) =K(s(z)—s)H(s(z)—s)G(s(z), s, s(z), s), where H is 
the usual Heaviside step function. We insert the ansatz into eq. ( |6.27f) and obtain 

= H(s - s)R{s - s) \d s d- s - 2 N9 1 G + 6{s - s)5{s - f)[G - 1] 

I (s + s) 2 ) 

+H(s - s)8(s - s)d s G + H(s - s)8(s - s)d- s G . (6.37) 

In order to solve this equation we have to find one solution G of the homogeneous 
equation (|6.3|) with the following properties: G(s — s, s — I) = 1, 9gG| s= j = and 
9 S G| 5= | = 0. It can easily be verified that a solution is given by 

e = ih { 2s{z) ' 8 + 8 - HijTW)^ + " 4 ' (6 ' 38) 

Thus we have determined an explicit expression for the Green function G = H(z— z)H(z— 
l)G. 



7 Related extremals 



Having determined solutions of the HJE ( j4.1|) combined with the IC (4.2) associated with 



a given extremal in terms of power series we now want to indicate how new extremals 
can be generated from a given one. 

In order to connect the functions Xn from equation (|4.20| ), with a one parameter family 
of extremals <p(z, z, u) embedded, we expand y = ((p—ipQ(z 1 z)) in the parameter u of the 
solutions of the equation of motion. Therefore we have to consider the two dimensional 
submanifold E:={(z, z, tp(z, z))} of the extended configuration space M.2+1 — {( z , z, <£>)}. 
The starting points of this calculation are the slope functions Q4.4|) : 



00 1 00 1 

d z (p{ Zl z) = d v s\ ¥= z = + E j-^-M > 9 *y = E -A+iy i+1 > ( 7 -i) 

i=2 I 2 i J- 8=1 % - 

00 1 00 1 

d-^(z,z) = d v S\^ = A 1 + J2j—TT;A l y\ => = £ 1 ^W +1 (7.2) 

i=2 i J- 8=1 

with y = (p — (p Q , A\ = dz^po and = <9 2 <^o- Expanding y in a power series of u 

00 _ u k 

y(z, z, u) = tp(z, z, u) - <p {z, z) =m^ A k (z, z)— (7.3) 

fe=0 K - 



yields: 



00 1 00 1 I 00 1 

^EV 1 = E 4i EV| , (7-4) 
fc=o i=i [j=oJ- 

00 1 00-1 [ 00 1 

^E> +1 = E 4i E^y^ . (7.5) 

fc=o ^ i=i I i=o J- 



Comparing the coefficients in Z + l-th order leads to: 

d z A t = AaA, + Inh(A a ,... J i4i +1 ;A ,...,Ai_ 1 ) 
d z Ai = A a Ai + Inh(Aa,...,Ai + i;Ao,...,A I _i) 



(7.6) 
(7.7) 



with the inhomogeneities Inh which depends on A4, Ai calculated above and the functions 
of lower order Aj,j < I. Making use of A 2 = dz\n(9) and A 2 = d z \n(9) we get: 



A; = el J -hih(A 2 ,...,A l+1 ;A ,...,Ai-i)dz + const 

= |y*-Iiih(A 2 ,...,i4i + i;Ao,... J A,_ 1 )dl + const} . (7.8) 



So the solutions of the Euler-Lagrange equations <p = y?o+I]fcLo AkU k+1 /k\ can be obtained 
by successive integration of the coefficients A;. Regarding the lowest order (1 = 0) d z A = 
A 2 A , d z A = A 2 A leads to: 

A = c 9, c EH, y = uc 9 + .... (7.9) 

Since 9 obeys the linear PDE (|4.12|) cq can always be absorbed into it. 9 is discussed in 
the ch. Applications above see ( |6.7|) , ( |6.8| ), ( |6.16|) and ( |6.17| ). 

Thus the one parametric family of extremals in the vicinity (u <C 1) of the 



original solution of the equation of motion tpo(z, z) is determined by: 



<f(z, z, u) = <p (z, z) + u9(z, z) 



(7.10) 



We compare our considerations with an expansion of the field variable (p within the 
Euler-Lagrange equation, (p is expanded in the parameter of the families of the extremals 
denoted as u: 



(p(z, z, u) = <p Q (z, z) + y(z, z, u) = <p (z, z) + u J2 ~r A " 



z,z)u 



n=0 



nl 



(7.11) 



<Po(z, z) is an arbitrary extremal. The potential is assumed to be a analytic function of 

(p: 



00 1 00-1 

v @) = E —wymote - v^r = e —d™vmw- ,„• 



m=0 



(7.12) 

Similar to the results of chapter 3, where we study the Hamilton- Jacob i theory, we 
obtain by comparing the coefficients of equal order of u n : 



u° : d z d z ip (z,z) + -d$V '(0)\^,=^ = 0, 



(7.13) 



which is fulfilled due to assumption and: 



u 



u 



n+l 



d z d- z A + -9?V(^ )Ao = 



(7.14) 



d z d- z A n + -5|y(^ )A n = InKAo.-A^, dlV(^)...d n ^ +1 V(^)) . (7.15) 
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The first one of these PDEs yields: 



(7.16) 



Because A has to fulfil the same linear equation as 9, which we introduced in the 
Hamilton- Jacobi theory ( |4.12|) , we are able to identify A and 9 with each other. Thus 
our result (|7.10|) is equivalent to a second variation of the action functional, which is 
commonly employed e.g. with semiclassical considerations 
tions 



13 1 and stability investiga- 
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The PDE ( [7.15| ) is analogous to the equation fl4.20p , which we obtained in the Hamilton- 
Jacobi framework. Both can be used to determine the fluctuations in a neighbourhood 
of a given extremal (po in every order of u. 

If one is able to find the general integral of the Hamilton-Jacobi equation and the inte- 
grability condition, the general solution of the Euler-Lagrange equation can be obtained. 



7.1 An example 

Here we would like to calculate the Hamilton-Jacobi functions S(z, z, <p), S(z, z, if) and 
a related family of extremals <f(z,z,u) according to the formalism developed in the 
chapters 4, 6 and 7. For this we choose the Sine-Gordon model with the 1-kink-solution 
<p = 4 arctan(exp(,z + z)). Though the functions S, S and the related extremals are 
determined perturbatively the corresponding formal series ( f4.5p , ( [7.3| ) can be obtained 
explicitly. 



7.1.1 The Hamilton-Jacobi potentials S, S 



The coefficients A n (z, z) and A n (z, z) ( |4.5| ) are determined by calculating the functions 
Xn(z, z) and x„(z, z) as solutions of the PDEs ( f4.20|) in every order n — 0, 1, 2, 



The coefficients A (z,z), A (z,z) can be calculated from ( |4.7|) , (|4.8| ) where the La- 
grangian C is given by 



i 



on the single extremal y?o- Then we obtain 

A {z,z) = -Acos(^y^ + d zXo(z,z)= A tanh (z + z) + d 2 Xo(z,z), (7.18) 
A (z,z) = -4cos(^) -d z xo(z,z) =4tanh(z + z) -d z xo(z,z). (7.19) 



2 

The coefficients of first order are determined by the embedding conditions: 

2 - 2 

A ± (z, z) = d 2 (p = —^j , Ai(z,z) = d z tp = —^j , (7.20) 

with the substitution l = z+z. Obviously this expressions for the coefficients Aq, Aq, Ai, 
At obtained are the general solutions of these equations (|7.18| ), (|7.19|) , (|7.20|) in contrast 
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to the following coefficients of higher orders of 

For the coefficients of the second order A 2 (z, z), A 2 (z, z) we choose the translation mode 
= l/cosh(Z) of Q 

A 2 (z, z) = d- z ln(0) = - tanh(Z) , A 2 {z, z) = d z ln(0) = - tanh(Z) . (7.21) 

The inhomogeneity of the equation for X3 (see ( ff.l7| ) and ( |4.22| )) always vanishes. There- 
fore, according to (|4.14|) , (|4.17j) and the discussion thereby we can choose Xs{ z iZ) = 
without any loss of generality. Hence the coefficients A 3 and A 3 are only determined by 
the function X3(z,z), which fulfill the inhomogeneous PDE (|4.20|) : 



tanh(Z) 



d z d- zX 3 - {2tanh 2 (* + z)- 1} X 3 = -^V{ V )\^ Q 6 2 = 2-—^ 
A special solution is given by X3 = ~ tanh(/)/(2 cosh(Z)), which yields: 



AJz,z) 



f-* It 



i 



2 cosh(Z) 







d 



A'3 



2 cosh(Z) 



(7.22) 



(7.23) 



For n = 4 the inhomogeneity of the wave equation (4.14) for the coefficient Xi ( |4.17D 
vanishes. Thus without any loss of generality we may set X4 = 0. The inhomogeneity of 
the equation ( [4. 14] ) for the function xa a ls° vanishes. So we may choose the translation 
mode again: 



Aa(z,z) 



^ ( ^ 



2 



e 



- tanh(Z) , A 4 (z, z) = ^d z 



- tanh(Z) 



(7.24) 



Regarding the series of coefficients we see, that those for odd and even indices contain 
the same functions 1/ cosh(Z) and tanh(Z), respectively. So we assume the same for the 
higher orders: 



.4 



2/i 



.4 



•l) n 



2n 



)2n-4 



tanh(Z) and A 2n +i = A 2n +i 



(-l) n 1 

2" cosh(/) ' «- 0,1, 2, 3, 



Inserting these expressions into the formal expansion ( f4.5|) we obtain: 



n n ,2n s 



S(z,z,(p) 



S(z,z,y) 



4 £ 

\n=0 

4 cos 
—4 cos 

(oo 
E 
71=0 



■l) n y 



n\2 2n , 
tanh(Z) + 4 sin 

y + y?o 



tanh(Z) +4 ^ 

\n=0 
1 



n„,2n+l ' 



11 



l 2 2n+l J cosh ^) 



2 J cosh(Z) 



2 

-i)V n 

n\2 2n 



4 cos 



—4 cos 



tanh(Z) + 4 sin 



+ d z Xo(z, z) = -4 cos ( - ) + d z xo(z, z) 



tanh(Z) + 4 iJ2 

\n=0 
1 



+ d z xo(z,z) 
.2, 

-i)V n+r 



n 



l 2 2n+l J C0S h(Z) 



2) cosh(Z) 



- d z xo(z,z) 



- d z xo(z, z) = -4 cos ( - ) - <9jXo(^, 2) 



(7.25) 

(7.26) 

0*Xo(M) 

(7.27) 



These solutions satisfy the HJE ( |4.1| ) and the IC ( |4.2|) . The embedding condition ( [4.4|) 
is fulfilled on the single extremal ipo by construction. 
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7.1.2 Embedded extremals 



The embedded extremals (p(z,z,u) can be determined by a straightforward integration 
of the eqs.: 

d z Lp(z,z,u) = $ = d ip S(z,z, ip = (p{z,z, u)) = 2 sin \^\ , (7.28) 
d z fi(z,z,u) = $ = d lf S(z, z, cp — <p(z, z,u)) — 2 sin (^j (7.29) 



which leads to 

dw 



2 sin (w/2) 



In arctan ( — 



(!)=*+ /® > ( 7 - 3 °) 



lnarctan(^) = * + /(*). (7.31) 



2 sin (w/2) V4, 

Obviously this system of algebraic equations can only be satisfied by the functions f(z) = 
z+ wand f(z) = z+ u, u = const, which gives the family of extremals: 

0(z, z,u) = 4 arctan (exp(z + z + u)) , (7.32) 

parametrized by one parameter u. 

However, here we would like to show how the embedded extremals can be calculated by 
the recursive formalism developed above, which is necessary if we are not able to deter- 
mine the DeDonder & Weyl Hamilton-Jacobi functions or the corresponding familily of 
extremals explicitly. 

First we calculate the first three orders (A ,A 1) A 2 ) of the expansion (|7.3|) of y(z,z,u) 
in u which will turn out to be sufficient to guess the general result for the embedded 
extremals y. From the eqs. (|7. 1|) , (|7.2j ) and (fPp we get 



d z y = A 2 y + ^A 3 y 2 + ^A 4 f + 0(u 4 ), (7.33) 
d- z y = Aiy + ^Asf + ^f + Oiu 4 ), (7.34) 



in which the formal expansion: 

00 _ u k 

y(z, z, u) = <p(z, z, u) - (p (z, z) = uY^ A k (z, z)— (7.35) 

k=0 

has to be inserted. For the first orders we obtain 

u l : d z A = )-d z (e)A , (7.36) 
u 

d- z A = ^d z (d)A , (7.37) 
.2 . a a _ 1 a ta\ a , 1 a ( *3\ a2 



^Ai = - Q dM^ + 2j£0. [fj A ° ' (7 ' 38) 
9mM = i^WAi + ^d- z Al , (7.39) 
.3 . a a _ o 1 a /aw o 1 3^*3,^3^ 



d z A2 = 2—d z {6)A 2 + 2—d z [y) A o + ^ ) A A, , (7. 10) 



9 Z -A 2 = 2—d z (6)A 2 + 2—d- z + ^ [f ) AnA, . (7. II ) 
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The solutions of these eqs. are determined up to a constant which may be absorbed by 
a redefinition of the parameter u. 



2 , tanh(Z) d / 1 \ . 

a 2 = ^ + 4 |_i^_^_^ mh2(() _i cosh2(() 



d 



2 



d/ 2 \3cosh(Z) # 



(7.43) 



From these coefficients we can can already guess the form of the coefficients of the higher 
orders: fc 

Ak = I+idi^ (c^sh(o) (7 ' 44) 

from which we obtain the embedded extremals: 

~ u k d k ( 2 \ [i™ u k+1 d fc+1 / 2 \ 

= ^ + f ^ dr - / ^Jm dl ' = 4 arctan (exp(/ + u)) - (7 - 45) 

Obviously we get a one parametric family of extremals satisfying the E.L. -equations, as 
well as the equations (7.28) and ( |7.29| ). It covers a strip of the extended configuration 



space: < ip < 2tt, z, z G 1Z. By translations ip^ip + 2n the whole 7Z 3 parametrized 
by z, z, (f G 7Z can be covered by families of these extremals with the exception of the 
parallel planes <f = 2kn, k = 0, ±1, ±2, . . ., which are solutions of the equations of motion, 
too. These are the so-called "vacuum" solutions in the Sine-Gordon theory. So this set 
of families of extremals, counted by the integer number k can be completed by these 
planes {p = 2kn, k = 0, ±1, ±2, . . ., so that the whole space 1Z 3 is covered by extremals. 



8 Wave Fronts 



In order to determine the wave fronts we have to turn to Caratheodory's framework, i.e. 
it is necessary to transform the DeDonder & Weyl Hamilton-Jacobi functions S(z, z, ip) 
and S(z, z, ip) appearing in the expansion ( ^4.5| ) and by the series ( |4.22| ) into those of 
Caratheodory S z (z, z, ip), S z (z,z,ip), namely S z and S z of eq. (|3.11|) . 



8.1 Caratheodory's Hamilton— Jacobi functions 

As stated in chapter 3 the two Hamiltonian formulations are algebraically equivalent for 
fields with only one field component. Thus the Hamiltonian density 7i and the momenta 
p, p are the same in both formalisms. Therefore we can use the equality of the momenta 
in order to determine S z , S z from the functions S and S calculated above. 
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We therefore return to the transversality conditions in ( |3.9|) , ( |3.11|) and obtain: 



OO 1 

d g S%S* -d z S z d v S' z = p = d v S = J2 1 A+iy l 

oo 1 

d z S'dyS* - d z S%S z = p=dyS = , £^A i+1 y i 



i=0 



U) 
5.2) 



Inserting the functions S" 2 and S z expanded in powers of y = (p — (p (z,z), like S and S, 



5.3) 



i=0 



i=0 



and comparing powers of y n yields: 



"-i-i ~~ E 

i=0 
n 



A 

An+l + E 



i=0 



n 
i 

n 
i 



d z A z - A* +1 d z(f0 
d- z A z - A z +1 d,<p 



A 



n—i+l 



d z A z - A z +1 d zVo 
d- z A\ - A z +l d m 



0. 
0. 



These equations determine the coefficients A^ +1 (z,z), A^ +1 (z,z) recursively: 



ip z A z )A z n+l — (d z Al)A z n+l — A n+ i + Inh" (Aq, Aq, . . . , A z n , A^) , 
{dzAl)A z n+1 — (d z Al)A z n+l = A n+ i + Inh^ (Aq, Aq, . . . , A z n , A%) ■ 



5.4) 
5.5) 



The functions Inh™ +1 and Inhg" 1 " 1 depend only on the coefficients of lower powers of y, 
namely Aq, Aq, . . . , A z n , A z n . They vanish for n = 0: Inh} = Inh^ = 0. For n = 1 we get 



Inh? = A\d z A\ - A\d z A\ and Inh^ = -A\d z A\ + A\d- Z A\ 
and for arbitrary orders n > 1: 



5.6) 



ra+l 



InhJ 



Inh™ +1 



n 



n-l 

-E 

i=i 

+Af cU* - A\d z A z n , 

n-l 



{ [9,A 2 - 4+i^o 



.4 



n— i+1 



d z A z - A z +l d z ^ 



n 
i 



+ E 

i=l 

-A\d z A z n + A\d z A\ 



{ [9 2 A| - Af +1 9^ 



.4 



n— i+1 



4i-i+l} 

(8.7) 



The determinant 



A = (d- z Al){d z Al) - (d z A z )(d z A z ) = C\ 



(p=(fi 



£o = d z (p Q dz<p - V((p ) 



•9) 



of this linear system of algebraic equations is assumed to be not zero. We consider 
only those extremals <po and regions in the Minkowski space where the Lagrange density 



(p=tp 



does not vanish. It makes sense to exclude those focal points and caustics, where 
v=!fi0 = 0, because the transversality relations of the wave fronts and extremals are 
violated otherwise Iffll. 
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To obtain the functions Aq and Aq we have to consider the zeroth order of the Hamilton- 



Jacobi equation in ( p. 11] ): 

(d- z Al)(d z A z ) - (d z A z )(d- z A z ) = d zVo d m - V(tp ) 



C . 



.10) 



If £o7^0 then one of the functions Aq, A$ can be chosen arbitrarily and the other one has 
to be calculated according to this PDE. For example we may choose Aq = z, A^ = J £ dz. 
In the case of the kink solution ipo = ±4 arctan(exp(2-|-,z)) for the Sine-Gordon theory 
we get AQ = 8tanh(z+z). 

The Hamilton- Jacobi equation of Caratheodory is satisfied automatically in any order 
of y n , n > 1, because the corresponding equation of DeDonder & Weyl is fulfilled: 

{d z s z d z s 2 - d z s%s z + n} = d zP + a- z p + d^n = d v {d z s + d- z s + n} = o . (8.11) 

The integrability condition holds too, because it imposes a constraint on the slope func- 
tions which are independent of the special Hamiltonian description for one component 
field theories. 



The linear system of equations for A z (z, z), A z (z, z) can be solved easily: 

-n+l 



A 



n+l 



A 



1 

1 



n+l 



(d z A z )A n+1 + (d- z A z )A n+1 \ + In]^ (A z , A z , 

- — 1 - n+l 

(d z A z )A n+1 + (d- z A z )A n+1 +Inh 2 (A Z ,A Z , 



(8.13) 



n+l n+l 

The functions Inh x , Inh 2 are linear combinations of the inhomogeneities ( |8.7| ) and 



^8j) may easily be determined. So we get e.g. for the coefficients A\ , A\ and A%, A?, the 
expressions: 

1 r,_ ... . ... [(d z A z )d z <p + (d z A z )d z <po} 



A\ 
A\ 



Co 
1 

To 



(d z J^)A 1 + (d z A z )A 1 
(d z A z )A 1 + (d z A z )A 1 



d z (p d z (po - V(ip ) 
[(d z A z )d- z ipo + (d s A*)d e <p ] 



and 



Ai 



.1 



1 

To 
1 

To 



(d z A z ) {A 2 -A\d- Z A\+A\d- Z A\) + (d 2 A z ) (A 2 + A\d z A\- A\d z A\ 
[d z A%) {A 2 -A\d- Z A\+A\d- Z AX) + (d- z A z ) (A 2 + A\d z A\- A\d z A\ 



(8.14) 
(8.15) 

(8.16) 
(8.17) 



The two coefficients A^ and A^ have to fulfil only one equation ( 8.10|) , because the 
closed 2-form J7 = dS' 2 A dS z is invariant under transformations S z — > S Z (S Z , S z ), S z — > 
S Z (S Z , S z ) with a Jacobi determinant that equals one. A family of wave fronts given by 



S z (z,z,Lp) 
S z (z 1 z,(f) 



£ -y A l (*> z ) (v 2 - <Po( z , Z )Y = O = COIlst - 
i=0 l - 

00 

J2 -i A i( z ' z ) (V 2 - ^oO, Z )Y = ° = const- 



l.tt 



(8.19) 



i=0 



is not changed by this transformation, but only reparametrized [^UJ a — > a = S z (a, a) 
a — > a = S z (a, a). 
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8.2 An explicit representation for the wave fronts 



In order to obtain an explicit expression for the one dimensional wave fronts z(a, a, if) 
and z(o~, a, (p) — cr, a fixed — we have to invert the relations ( |8.19| ) which can be regarded 
as the defining equations for these functions z(<p), z{ip). Here they are assumed to be 
analytic functions of ip: 

DO j OO 

z = J2- ai (a,a) (y = (p-(p (a,a)y , z = J2 °) (v = <p-tPo(?, °)Y > ( 8 - 20 ) 



i=0 



11 



i=0 



II 



i.e. they can be expanded in powers of the difference y — (p—<po(a,a), where 0o(a,a) = 
ipo (z(a, a), z(a, cr)) denotes the extremal in terms of the variables a, a. Obviously there 
is a difference between the quantities y = (p—{p (z, z) and y = <p—(p (a, a) . Thus in order to 
determine the coefficients a n , a n we have to insert this series into the defining equations 
for the wave fronts ( 8.1 8|) and (|8.19 ): 



oo i / oo yfi oo y^ra \ 

o- = -T A i\ z = a o(c, a) — r"n(c r , a), z = a (cr, a) —,a m {o-, a) ,(8. 21) 



i=0 



n=l 



nl 



m=l 



ml 



a 



oo i / oo 



i=0 \ n=l 

where the variables y have to be expanded in powers of y, too: 



nl 



oo ~m \ 

a (a,a)+Y^ —,a m {a,a) \ (8.22) 



y = (p-(p i z = a (a,a : 



+E 

n=l 



nl 



a n {a,cr), z 



m=l 



a (a,a) +V —a m (a,a) . (8.23) 
ml ' 



The wave fronts on the extremals ao(a, cr), cxq(o~, cr) are determined by the zeroth order 
of the eqs. (|OTD , QOg ) and ( g23[) : 



a = (2 = a (o", a), z = a (cr, cr)) and cr = [z = a (o~, a),z = a (cr, cr)) . (8.24) 

For the Sine-Gordon theory we obtain from the choice we have made for the functions 
Aq, Aq\ z = a (cr, cr) =cr, z = a (o~, cr) =artanh(cr) — cr. 

Locally the functions «o(cr, cr) and ao(a,a) are determined uniquely, because the Jacobi 
determinant A = {dzAfy(d z Afy— {d z Al){d z AQ) =C\ v==l p of this transformation z,z — > a, a 
does not vanish, as assumed. 



Expanding the expressions ( p. 21 ) in y n , inserting the series ( |3.23| ) and separating the 
coefficients a n (a,a), a n (a,a) we obtain a pair of linear algebraic equations for these 



coefficients, that can be solved like the system ( |3.4|) , (|8.5| ). After inserting the coefficients 
Ai, Ax, . . . , A n , A n of the DeDonder and Weyl Hamilton-Jacobi theory we get the final 
result for n > 2: 



f \n— 1 



a, 



a, 



n_in- 

'-Cor- 1 

m +1 



<A)) 2 ^n) +Inh(y4i, Ai, . . . , A n -i, A n _i, deriv.), (8.25) 



(v(ip )A n -(d z ip ) 2 A n ) +Inh(Ai, Ax, A n _i, A n _ l5 deriv.). (8.26) 
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The term "deriv." in the inhomogeneities denotes the derivatives of A z , Af, %— 1, . . . , n — 
1 with respect to z and z. The functions A z nl A z n and their derivatives are short-cuts of 
e.g. d z A z n = d z A z n \ z=ao ^)^=a {(T,a) and thus depending on a, a only. 
Obviously a not degenerated transformation z,z — > a = S z ,a = S z only exists if the 
Lagrangian £ an d the Hamiltonian densities Ho on the extremals do not vanish. 

The coefficients of zeroth order are given in ( |8.24j ), those of the next two powers in y 
are: 

/ — \ d z ip dgtpo _ d z cp d z <p 

ai{a,a) = — — = — — — r , a\{a,a 



and 



Ho d z (p d z (po + V((fo) ' H d z (p o d z (p + V(ip ) ' 

(8.27) 

a 2 {a, a) = V [-C d s ln(6») + d- z (d z ip Q d E ip )} 
rio 

d^fo^fo ln {9) + ^ (dz(podg(po)] + ^o^v Wj (8 .28) 

n rto 

a 2 (a, a) = V ^ [-£ d z ln(0) + d z (d z (p d z (p )} 

- dz(f ^f° [-C d z \n(6) + d- z {d zVQ dm)\ + • (8-29) 

rto rto 

Thus the wave fronts z(a, a, (p) and z(a, a, (p) can be calculated from the coefficients A n , 
A n of DeDonder and Weyl's Hamilton-Jacobi framework and the coefficients Aq, Aq of 
Caratheodory's one. Caratheodory's coefficients A z n and A x n , n > 1 are not necessarily 
be involved in the final expressions ( |8.25| ), ( |8.26| ), since they can be substituted by those 
of DeDonder and Weyl — determined by Backhand transformations. 



8.3 An alternative way to determine the wave fronts 

If one is not interested to get an explicit representation of Caratheodory's Hamilton- 
Jacobi functions S z (z,z,(p), S z (z,z,{p) it is possible to get the coefficients a n , a n of 
the last paragraph much more easily. However, this method can only be applied, if the 
Hamiltonian density H does not vanish! If H = we have to calculate Caratheodory's 
Hamilton-Jacobi functions S z (z,z,ip), S z (z,z,ip) as described above. 

For nonvanishing Hamiltonian densities H we may change the independent variables 
z, z to a — S z (z,z,ip) and a = S z (z,z,<p), while the field ip remains unchanged. The 
functional determinant of this transformation is just H according to Caratheodory's 
Hamilton-Jacobi equation ( |3.11| ). Then the wavefronts z = z(<p, a, a), z = z(<p,a,a) can 



be determined explicitly from the equations 

d v z{(p, a, a) = , d v z(ip, a,a) = — . (8.30) 

These equations are obtained by comparing the coefficients of the wedge products dip Ada 
and dip A da in the equation (|3.11|) , if the variables a = S z and a = S z and ip are regarded 
as the independent ones P2J. So we can e.g. immediately determine the coefficients of 
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the first order ai and a\ ( |8.27D . The remaining coefficients a n and ot n in (|8.26|) and 
( p. 25 ) can be determined by expanding the DeDonder and Weyl momenta p, p and the 



Hamiltonian density H in powers of y. 



8.4 An example: the Sine— Gordon theory 

We return to our example in chapter 7 in order to illustrate the formalism discussed 
above. 



8.4.1 Caratheodory's Hamilton— J acobi functions 

To obtain Caratheodory's Hamilton-J acobi functions S z (z,z,tp), S z (z,z,f) from those 
of DeDonder & Weyl, ( |7.26|) and (|7.27|) , we first have to solve the equation (|8.10|) in 



order to get the coefficients Aq(z, z), Aq(z, z) of order y°. Because in our example the 
Lagrangian density depends only on l = z + z on the extremals it is useful to transform 
the independent variables z, z into l = z + z and l = z — z 

(djA^Al) - {djAimAD = l -C Q = — * . (8.31) 

Due to the invariance of the 2-form f2 = dS^ A dS z with respect to the transformations 
S z — > S Z (S Z , S z ), S z — > S Z (S Z , S z ) with a Jacobi determinant equal to one we can choose 
any single solution of this equation without any loss of generality, e.g. 

A z = A Z {1) = 4 tanh(z + z) , A\ = A\ = -I = z - z . (8.32) 



Inserting this result into the equations ( |8.14| ) and flS.15| ) we get the coefficients of first 
order 

cosh(/J 



which leads, using the eqs. (|8.12|) , ( ^.13j ), (|3.16|) and ( |3.17|) , to the coefficients of the 
second and the third order in y 

^ = ^(/) = -tanh(/), Al = A z S) = -^-- , A* = A* = 0. (8.34) 

Inspecting the series of these coefficients we get an ansatz for those of all orders y n , 
n > 1 

^n = ^L(0=4^tanh(Z), Aln _ x = AL _ l{ i) = 4 Hj" * , ^ = 0, 

(8.35) 

from which we obtain Caratheodory's Hamilton- Jacobi functions S z (z, z, </?), S z (z, z, tp) 
by determining the sum ( p.3|) : 

S(z,z,<p) = 4 ( V ^y n ) tanh(Z) + 4 ( V ^^- ) —^ttt , 

= 4 cos (|) tanh(/) + 4 sin (f ) = -4 cos (f ) , (8.36) 

S 2 (z,z,(p) = z-z. (8.37) 
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These functions satisfy Caratheodory's Hamilton-Jacobi equation ( f4.1|) and the integra- 
bility condition ([4.2|) as well as the embedding conditions (fOf), ( |4.4|) . The solutions 
( p.36| ), ( |8.37| ) for S z , S z are the only ones provided that we have fixed the coefficients 
Aq and Aq, because the eqs. (|8.12| ), ( |3.13|) , ( fB.16| ) and (|3.17|) lead to unique solutions . 



The general solution of the system of PDEs ( (4.2| ) and ( (4.4| ) is obtained by apply- 

ing arbitrary smooth transformations S z — > T,(S Z , S z ), S z — > T>(S Z , S z ) to the functions 
( gg§ ), ( gg7| ) with a Jacobi determinant <9i£<9 2 £-<9 2 £<9i£ = 1 

S^z, 2, </?) = £ (2 — z, — 4cos (<f/2)) and S*(z, z, = £ (z — z, — 4 cos (<f/2)) . 

(8.38) 



8.4.2 An explicit representation of the wave fronts — the singular case 



Caratheodory's Hamilton-Jacobi functions S z , S z are given now by the equations (|3.38|) 



If we would like to have them in the explicit form z(a, a, if), z[o, a, ip), we have to invert 
the equations (|8.18|) , (|8.19|) . Doing this we get only the zeroth order coefficients 



ao(a, cr) = - ^artanh— + tjj and ao(a, a) = - ^artanh— — o^j (8.39) 

without any difficulties. All the other coefficients a n , a n , n > 1 do not exist, because 
the Hamiltonian density TC = pop + V (if) with V = 2(cos(^) — 1) vanishes on the given 
extremals <p (see (ggg ) and (|Q6D )! 



From the results (|8.36|) , ( p. 37] ) it becomes obvious why this happens. The function ( |3.37|) 



depends only on the difference z—z whereas ( ^.36 ) does not depend on the variables 



z, z 



at all! Therefore it is impossible to invert the equations ( |3.18|) , (|8.19|) to get the functions 
z[o, a , f) and z[a, o, if) at every point (cr, a) of the parameter space T = {(a, a)}. Hence 
the wave fronts are only defined at points (a, a, if) for which 

<p = g(o-) = 2 arccos (a/A) , -4 < a < 4 (8.40) 

holds. The real parameter a remains arbitrary. Nevertheless the wave fronts are one 
dimensional straight lines parallel to each other: z — z = o = const.. They cover the 
extended configuration space M.2+1 = {z, z, ip} as required, but the wave fronts cannot 
be given as functions z[a, a, ip) and z[a, a, ip) because they are parallel to the surfaces 
if = const.. 

Due to the transversality condition (|3.18|) the 2-parametric family of 1-d wave fronts 
(see Fig. 1) 

z — z = o = const. , ip = g(<j) = const. (8-41) 

and the 1-parametric set of the 2-d extremals (p = 4 arctan (exp(/ + u)) intersect each 
other transversally everywhere in M2+1 = {z,z,ip}, because the Lagrangian density C 
differs from zero in M.2+1 contrary to the Hamiltonian density which vanishes everywhere 
in M.2+1- The angle l{w,e) between the basis vector W of the tangent space TpW of 
the wave front 

w = 2 sin(yV2)(d z + 3?) (8.42) 
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and an arbitrary nonvanishing vector e in the tangent space T-pE of the extremal para- 
metrized by the variables A, A 



e = Xd z + Xd B + 2(A + A) sm{ip/2)d v , with |A| + |A| > (8.43) 
at the point P — (z, z, <p) is given by 

f (w,e)\ ( A + A \ . A , 

L{w,e) = arccos — rr—r = arccos _ _ = , (8.44) 

V y/2 J\ 2 + A 2 + 4(A + A) 2 sin%/2) j 



v \ w \ \ e \ 

the minimum of which with respect to a variation of the parameters A, A gives the angle 
under which the wavefront intersects the extremal at P: 

Z(W, £) = rninZ(u>, e) = arccos | , | = arccos | = | , (8.45) 

^l + 8sin>/2)J Vv^V2/ 

Z(W, £) is constant on a given wave front a = const., a = const.. This angle Z(W, £) 
takes its minimum Z(W, S) = for the values (p = (p^ = 2kir, k = 0, ±1, ±2, . . ., namely 
Z(W, S) = 0. So in the limit z + z = I — > ±oo the transversality condition is vi- 
olated, since the wave fronts and extremals ( |7.45| ) are asymptotically parallel at all 



points P in the planes {z, z G 7Z, (p = 2kn] e M.2+1 parametrized by the integer num- 
ber k = 0, ±1, ±2, . . .. Any of these planes separates the two one-parameter families 
<^ = arctan(exp(/-Ki))+2(/c— l)7r and <^ = arctan(exp(/-|-M))-|-2A;7r. Every extremal of these 
families "touches" the plane (p = 2kir in the limits I — > —00 or / — > 00. The maximum 
value for Z(W, £) = arccos(l/3) indicates that the wave fronts and the extremals are 
never perpendicular to each other. 

Notice, that the singular situation discussed above is essentially a coordinate singularity 
resulting from the singular change of variables z, z to a = S z (z, z, ip) and a = S z (z, z, ip). 
It shows that we cannot choose the parameters a, a, tp as independent variables to rep- 
resent the wave fronts. The wave fronts are not degenerated at all. 

However the transversality between the wave fronts and extremals is violated in this 
case on the boundaries of the regions A4k = {2kn < ip < 2(k + l)n, z,zE 71}, where the 
families of extremals ( |7.45|) are defined — a singularity that cannot be circumvented by 
a coordinate transformation. 

These results can be obtained in a straightforward manner by using the DeDonder-Weyl 
Hamilton-Jacobi functions S(z,z,<p) and S(z, z, <p) from ( [7.26| ) , ( |7.27| ) . 



The Hamiltonian density Ti. vanishes on the extremals <p(z, z, u) of ( |7.45|) 



U = d^Sd^S + V(<p) = (2 sin(<^/2)) 2 - 2(1 - cos(») = . (8.46) 
Notice that the "usual" canonical Hamiltonian density: 

n can .(x,t) = -7T 2 + -{d x ^f + (1 - cos(p)) with 7T = - (8.47) 

does not vanish on the kink solution ipo, because the Legendre-Transformation is applied 
only with respect to the time derivatives d t <p of the field ip and does not affect its spatial 
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derivative d x tp. This underlines again that the energy density H. can , is different from the 
"covariant" Hamiltonian density 7i Q3.5Q , that we use in this article. 



Caratheodory's Hamilton-Jacobi equation ( j3.11|) on the family of extremals ( |7.45| ): 

d z S z d- z S~ z - d z S z d- z S z = -H = , (8.48) 

the solutions of which determine the wave fronts transversal to the extremals we are 
interested in, shows that the change of the variables z, z to o = S z (z,z,ip) and a = 
S z (z, z, tp) is not a regular one, because 7i is nothing but the functional determinant of 
this transformation. 



The general solution of equation ( |8.48| ) is: 

S' = V{f(z,z),<p), S* = &(f(z,z),<p), (8.49) 

with an arbitrary smooth function f(z,z). Inserting this result into the conditions on 
the momenta ( |3.11|) : 

p = 2sin(vV2) = d z S z dvS z - d- z S z d v S z = +(dj) (d^d^ 2 - d v X z d f E*), (8.50) 
p = 2sin(^/2) = 8 z S z d v S' z - d z S z d v S z = ~{d-J) (dfE'dyE* - d v E'd f E s ) (8.51) 

shows, after subtracting these PDEs from each other, that the function / depends on 
the difference l = z— z only: f = f (z—z) . Due to the invariance of the basic 2-form ( |3.10| ) 
with respect to symplectic transformations we may choose / = z — z, S 2 = — 4cos(<£>/2) 
and S z = z— z without any loss of generality. This result coincides with (|8.36|) and (|B.37|). 



8.4.3 An explicit representation of the wave fronts — the regular case 

If we would like to circumvent the spurious singularity discussed above we have to ensure 
that the Hamiltonian density does not vanish in M.2+1, e.g. by adding a global constant 
Co to the Lagrangian density: 

£ =>. C = d z ^d- z y + 2(1 - cos(v?)) + c , (8.52) 

which does not influence the equation of motion. More generally, we can add any exact 
two form r = d(/(z, z)dz+f(z, z)dz) to the basic form (|3.2j) , Q — > Q+T, without affecting 
the momenta p,p, the slope functions v,v and the equation of motion, but modifying 
the Hamiltonian and the Lagrangian densities Tl, C: 

H - H + d-J-dJ, £ ^ £-d- z f + d z f. (8.53) 

We see that the form of the wave fronts is influenced by adding such a term T while the 
family of extremals ( |7.45| ) remains unchanged — a property which is known in mechanics, 
too. For the sake of simplicity we choose f = 0, J' = Cqz which shifts £ and Tl merely by 
a constant c . 

With the exception of the coefficients A , A of zeroth order, which are of no interest in 
determing the embedded extremals and the corresponding wave fronts, the DeDonder 
and Weyl Hamilton-Jacobi functions S(z, z, ip) and S(z, z, ip) and the resulting momenta 
p = d ip S, p = d !fi S are not affected by this shift contrary to the wave fronts. 
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Using Caratheodory's Hamilton- Jacobi equation ( 3.11]) : 

d z S z d- z S~ z - d z S~ z d z S z = -H = c (8.54) 
which represents nothing but the determinant of the linear system of equations (see 

(Q) 

p = 2 sin(^/2) = d- z S~ z d v S z - d- z S z d v S z , (8.55) 
p = 2 sin(p/2) = d z S z d v S~ z - d z S~ z d v S z . (8.56) 

For the functions d v S z and d v S z we get a system of decoupled linear PDEs for the 
functions S z , S z 

(Hd v + pd z + pd- z ) S z = {-c d tp + 2sm{ip/2)[d z + d z \) S z = , (8.57) 
(Hd v + pd z + pd z ) S z = (-c ^ + 2sim>/2) [d z + d z }) S~ z = 0. (8.58) 

The general solution: 

S z = S(-co2/2 + 2cos(^/2),-c z/2 + 2cos(^/2)) , (8.59) 
S z = S(-co2/2 + 2cos(^/2),-c z/2 + 2cos(^/2)) (8.60) 

of this system PDEs is obtained by the method of characteristics. The functions E, E 
have only to satisfy Caratheodory's Hamilton- Jacobi equation ( 3.11| ). Taking the in- 
variance with respect to symplectic transformations into account we may choose the 
functions 

S z = V2 [z - (4/c ) cos(p/2)] and S z = V2 [-c z/2 + 2 cos(<^/2)] (8.61) 

without any loss of generality. 
Obviously the transformation 

(z,z) ^ (<r = S'{z,z,<p), a = S- z (z,z,v)) (8.62) 

exists now, leading to an explicit representation of the wave fronts: 

\f2(j 4 a 4 

z(a, a, ip) = 1 cos(<^/2) , z(a, a, ip) = —= H cos(<^/2) (8.63) 

Co Co v2 c o 

with some properties different from those obtained in the "singular" case: the 1-d wave 
fronts are not straight lines in the extended configuration space. Subtracting the eqs. 
QB.63 ) from each other shows that they lie in planes parallel to the 1-axis like in the 
singular case. 

Similar to the "singular" case the angle l(w,e) between the basis vector of the tangent 
space 7pW of the wave front 

w = 2 sin(vV2)(c) z + d g ) - c d v (8.64) 

and an arbitrary nonvanishing vector in the tangent space T-pE of the extremal 

e = \d z + \d z + 2(A + A) wy[tp/2)d v , with |A| + |A| > (8.65) 
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at the point (z,z,ip) is given by 



f(w,e)\ ( 2(l-c )(A + A)sim>/2) 

Z(w, e) = arccos -j — — =arccos 



\w\\e\ 



4 + 8 sin%/2WA 2 + A 2 + 4(A + A) 2 sm 2 (cp/2) / 



the minimum of which with respect to a variation of the parameters A, A gives the angle 
Z(W, £) in which the 1-dimensional wavefront intersects the 2-dimensional extremal at 
P: 

, . ( 2v / 2|l-cbl|sin(p/2)| \ 
Z(W, S) = minZ(w, x) = arccos 1 yyi n . (8.66) 

V ^jcl + 8 sin 2 (yV2) VI + 8 sin 2 (<^/2) / 

Obviously in the case cq = 1 the wave fronts cross the extremals always perpendicularly, 
because Z(W, £)=ir/2 holds everywhere in .M2+1 (see Fig. 3). 

For Co 7^ 0, 1 the maximum value Z(W, £) = 7r/2 is taken only for ipj. = 2kir, k = 
0, ±1, ±2, . . ., contrary to the singular case cq = 0, where for these values of the field 
variable the minimum of the angle Z(W, £) vanishes, indicating that the transversality 
relation between the wave fronts and extremals is violated there. The angle Z(W, £) has 
a local minimum at those points where | cos(Z(W, S))\ becomes maximal: 

|cos(Z(W i g) mml )| = 2v ^ / |1 ~ C ° l , |cos(Z(W,£) mm2 )| = i^4. (8.67) 

3^/c 2 , + 8 1 + N 

We find the first minima Z(W, £) m i n i on planes ip = {2k + 1)tc, k = 0, ±1, ±2, . . . where 
the cosine of the field variable vanishes. They do exist for all values of the parameter 
CqEIZ \ {0}. For c 7^— 8 the angle l(W,£) m i n i differs from zero, i.e. the transversality 
relation is fulfilled. If Co = —8 both minima (|3.67|) coincide, therefore this case is discussed 
below. 

The second minima l(yV,£) m in2 exist only in the range |cq| < 8, Co 7^ 0. These minima 



lie at values of the field variable: ip = 2kTc ± 2 arcsin(y |co|/8), k = 0, ±1, ±2, .... If Co <0, 
this angle vanishes, whereas it differs from zero for all co>0. 

It results that in the case Co > or Co < —8 the wave fronts and extremals are never 
parallel to each other (Z(W, £) m mi,2 0)- This guarantees the transversality relations 
everywhere in the extended configuration space M.2+1 — even on the planes ifk = 2kii, 
k — 0, ±1, ±2, . . . that separate the one parameter families of extremals — contrary to 
the singular case. This result coincides with the fact that both the Lagrangian and the 
Hamiltonian densities do not vanish in M.2+1 for co>0 or co<— 8. 

For — 8 < c < the minimum Z(W, £) m in2 is equal to zero. This happens just for the 
points on the plane 



2A;7r±2arcsin(y / -c /8), k = 0, ±1, ±2, . . . (8.68) 



where the shifted Lagrangian density ( |3.52| ) vanishes, as expected from the transversality 



relation (|3.18|) (see Fig. 2). The case c = — 8 is a special one: here the two planes ( |3.68|) 



that exist in every region 2kir < ip < 2(k + z,z G 7Z, k = 0, ±1, ±2, . . . coincide. 
Hence we get only one plane f = 2kn + 7r in the range 2fc7r < ip < 2(k + 1)tt, where the 
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transversality relations are violated due to a vanishing shifted Lagrangian density. 

The dependence of the angle Z(W, £) on the constant c shows, that the geometrical 
properties of the wave fronts and even the transversality relation may be affected by 
changes of the Lagrangian density that do not influence the equation of motions and the 
corresponding extremals at all. 

Now we would like to show how these results can be obtained using the recursion formulas 
given in sections (8.1) and (8.2). We first calculate Caratheodory's Hamilton- Jacobi 
functions S x (z, z, <p), S z (z, z, ip) from those of DeDonder & Weyl, ( |7.26|) and (|7.27| ). For 
the sake of simplicity we choose the shift Co = —2 in order to eliminate the constant 2 in 
the Lagrangian of the Sine-Gordon model C — d z (pd s (p +2(1 — cos(<f)). 

The coefficients Aq(z,z) and Aq(z,z) can be chosen as special solutions: 



A* = A z (l) = 4tanh(Z) -/ 



A z 



I 



(8.69) 



of the equation flS.10|) written in terms of the variables l = z + z and l = z — z. 

1 . 



{d l Al){d l Al)-{d l Al){d l Al) = ^ ' 



2 cosh(Z) 2 

The equations (|8.14|) and QS.15j ) yields the coefficients of first order: 

2 



(8.70) 



A{ = A{{1) 



cosh(/) 



A\ = 



(8.71) 



and the relations ( |8.12j ), ( |8.13| ), ( |8.16| ) and ( |3.17|) lead to those of the second and the 
third order in y. They coincide with those of the singular case (|8.34j) since the ratios 



(OzAq)/ C and (dzA^/C are the same in the regular and in the singular case. The same 
holds for the coefficients A z n {z, z) and A z n {z, z) of higher orders in y. So we get for n>0: 



A 



2n 



2 2n 



■ tanh(/) , A 



2n-l 



^n-l(0=4- 



-l) n 



2 2 "" 1 cosh(Z) 



At = 



(8.72) 

which combine the expressions ( |3.69| ) yield the Hamilton- Jacobi functions S z (z,z,(p), 
S s (z,z,<p): 



S(z,z,cp) 



S z (z,z,(p) 



(oo 
n=0 



-i)V n ^ 

n\2 2n 



tanh(/) + 4 ^ 

\n=0 
1 



(-i) n y 



n„,2n+V 



4 cos - tanh(/) + 4 sin - 



2 J cosh(Z) 



n ] 2 2n+l 

-1 = - 



cosh(/) 



z — z . 



4cos( ^j-z-z, (8.73) 
(8.74) 



Both functions S z and S z satisfy Caratheodory's Hamilton-Jacobi equation ( |4.1| ), the 
integrability criterion ( [4.2D and the embedding conditions (|4.4| ). They are related to the 
solutions (|8.61| ) by the transformation 



S z 



S z 



~V2 



{S z + S 2 ) , s~ z -> 



V2 



(S~ z - S z ) . 



(8.75) 
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The wave fronts z(a, a, if) and z(a, o~, ip) can be determined recursively following the 
ideas discussed in the section (8.2): 

The coefficients ao(a, a), ao(a, a) of zeroth order in y are given by inverting the eqs. 



a = Aq [z = ao(a, a), z = ao(a, a)) = 4 tanh(a + ®o) — («o + ®o) , (8.76) 
a = Aq(z = a (a,a),z = a (a,a)) = -a + a , (8.77) 

which yields: 

a = ±(f- 1 (a)-a) and a = ~CTV) + *) , (8-78) 

where the symbol f~ l denotes the inverse of the function f(x) := tanh(x) — x. The 
coefficients of y, y 2 and y 3 are given by ( gg§ , ( |8~27| ), ( ggg ) and ( ^29|) : 

ft cosh(/ ^a)) ft cosh(/ ^cr)) 

« 2 (o") = ~ tanh(r 1 (o-)) , a 2 (a) = ~ tanh(rV)) , (8.80) 

"3(0") = -- r77Zi7 7x » "3(0") = -- tjjzu \\ ■ (8.81) 

4cosh(/ i (cr)) 4cosh(/ 1 (crj) 

This leads to a general ansatz for the coefficients a n and a n , n>l: 

(— l) n (— l) n 1 

a 2ri = a 2n = 2-^- tanh(/- 1 (o-)) , a 2n ^ = a 2n ^ = 2 — CQsh(/ _ 1((j)) > (8-82) 

which give the wave fronts 

/ 00 /_i\n 2n\ / 00 (_-i\n 2n+l\ i 

*(<^) = 2 tanh(fV)) + 2 fc ( ' " 1 



n\2^ J w v " V^o nl22n+1 J coshC/- 1 ^)) 

-I., , 1 , ,• I .. — \ . ..... (r\ 1 



-2tanh(/-») + - (/ -» - a) = -2 cos (JJ - -(a + a), (8.83) 
z(<7,<7 )¥ >) = -2 cos (I) -\{cr-ar). (8.84) 



which are related to the representation ( |8.63| ) by the transformation (|8.75|) applied to 
S z = a and S z = a, which reparametrizes the family of wave fronts only. 



9 Conclusions 

Within the manifest covariant Hamilton-Jacobi canonical frameworks of DeDonder & 
Weyl and of Caratheodory we have investigated relations between families of extremals 
and Hamilton-Jacobi wave fronts for 2-dimensional one component field theories. This 
is of interest, since the dynamics of fields can be described either by the Euler-Lagrange 
or the Hamilton-Jacobi equations supplemented by the integrability conditions. 

We developed a formalism to solve the DeDonder & Weyl Hamilton-Jacobi equation and 
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the integrability condition perturbatively by expanding the Hamilton- Jacobi functions 
in powers of the field variable. Starting from a single given extremal it is then possi- 
ble to calculate new ones from it by using the two DeDonder & Weyl Hamilton- Jacobi 
functions. 

This formalism is useful especially for investigating extremals in the neighbourhood of 
known extremals in several 2-dimensional field theories: the massless and massive Klein- 
Gordon, the Sine- and Sinh-Gordon, the Liouville as well as the 4 theory. In the Sine-, 
Sinh-Gordon and in the 4 theory we have studied the embedding of topologically non- 
trivial soliton solutions. This approach is related to the usual stability investigations of 
solitons where perturbations are considered which are a product of functions depending 
on the time and the space variables separately. We determined the general solutions of 
the equations of the second variation by using Backlund transformations. In this manner 
we have obtained all the extremals in the vicinity of a given one. 

Calculating proper Hamilton- Jacobi wave fronts makes the use of Caratheodory's Ha- 
milton-Jacobi functions S z (z, z, (p), S z (z, z, (p) necessary. Solving Caratheo's Hamilton- 
Jacobi equation is considerably simplified by using the corresponding DeDonder-Weyl 
Hamilton- Jacobi functions. One first obtains the wave fronts as equipotential surfaces 
S z = a = const., S z = a = const, in an implicit form. If the transformation of variables 
(z, z, ip) — > (a, a, ip) is a regular one, i.e. if Caratheodory's Hamilton density does not 
vanish on the family of extremals under consideration, we get an explicit representation 
z = z(ip, cr, a), z — z(<f, a, a) for the wave fronts. 

These general results have been applied in detail to a special single kink solution of the 
Sine-Gordon equation. After calculating the DeDonder & Weyl Hamilton-Jacobi poten- 
tials 5*, S we obtained a corresponding one-parameter family of embedded extremals: 
a field of kink solutions of constant energy covering the extended configuration space. 
From the functions S, S we have determined Caratheodory's Hamilton-Jacobi functions 
S z , S z explicitly. The wave fronts have been determined for the singular (7i = 0) as well 
as for the regular case (TCy^O). In addition the transversality conditions between the 
wave fronts and the embedded extremals have been analyzed. 



A Appendix 

Backlund transformations are employed to map the integral submanifolds of the PDE 
( |4.12| ) to those of the Klein-Gordon or wave equation, the general solutions of which are 



known. In this chapter we determine all the functions dfV((po{l)) with l = z+z that allow 
this map by one BT. But unfortunately it turns out that the nonconstant coefficients 
dfV((p (l)) in ( |6.12|) for the 4 models do not belong to this class. Thus, we need 

at least two BTs to connect these two equations to one PDE with constant coefficients. 
Moreover we have to discuss which relations of type ( }4.12|) can be reduced to free field 
equations by one or a finite number of BTs at all. 

We show that there is no BT which relate the PDEs (|6.11|) , (|6.12|) and d z dzQ = m 2 9, 



m 2 e 1Z. Even if we make the more general assumption that the transformed function 
9 obeys the eq. d z d 2 = m 2 (z,z)8 with a nonconstant but separable coefficient m 2 : 
m 2 = B(z)B(z) the equations ( |6.11| ) and ( |6.12j ) for the 4 theories cannot be reduced by 
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one BT to such an equation for 6 and thus not solved like the PDEs ( |6.3| ), (|6.5|) , ( |6.6|) 
for the Liouville, the Sine- and the Sinh-Gordon models. The reason is that equation 
d z d z 8 = m 2 (z,z)8 can be reduced to a Klein-Gordon one by a suitable transformation 
of variables z — > B(z) and z — > B{z) [18[]. The equations (|6.11|) and (|6.12|) can be 
transformed by two BTs to Klein-Gordon equations. This result can be generalized to 
PDEs of the following type: 

d z d- z 6 = {n[n + % 2 (7 = z + z) + a}6 n = l,2,..., aeTZ, (A.l) 

if 7] is a solution of dir] = br] 2 + c with 6 2 = ±1 and c G TZ by h BTs. Especially the PDEs 
( |5 . 1 1| ) , (|6.12|) can be obtained by choosing b= —1, c = 1 and n = 2. We start from eqs. 
( [5.12| ), assuming m 2 — > m 2 (Z) to be a function of l = z+z and rename d 2 V((po(l)) as 2v(Z). 
From the equations ( |5.12|) we infer 

^ = ^o + ^)+A(^), Ml) = lfdl'fdl"(m 2 (l")+v(l")), (A.2) 
^0 = ^o + (^o + 5 2 A)(9^o + ^)-m 2 (/) (A.3) 

with two arbitrary functions A and A. We are interested in the relation between m 2 (l) 
and v (I). First we have to determine the functions A and A, which depend only on z 
and z respectively. Differentiating the PDE (|A.3|) with respect to l = z — z leads to: 



d^o [d 2 z A - d£A] = diAd z A - d 2 z Ad- z A , (A.4) 

because ipo depends on / only. Now, if d 2 A = or d 2 A = 0, we obtain from eq. ( |A.4| ) 
that either d 2 A = and d 2 A = (case I) or diipo = —ao (case II). The case II leads to 
dfipo — and using relation ( |A.3| ) gives m 2 (l) = 0. The definition of tpo finally yields 
v = 0, which is of no interest for us. Without any loss of generality we have chosen: 
d 2 z A = Q => d z A = a , a E TZ. 

If = d 2 A-d 2 2 A (case III), then it follows that either d 2 A = d 2 s A = or d z A = d z A = a = 
const., the second case being a special case of the first one. So we have d z A = a\ and 
d z A = a\ with £ 1Z. Thus case III is contained in the first one. 

If d 2 A—d 2 A^0 (case IV) we are able to divide the relation (|A.3|) by it and to apply the 
operator dj once more: 



= {d 3 z Ad 2 z A + d 2 z AdlA)(d z A-d- z A) -2d 2 Ad 2 z A(d 2 A-d 2 z A), (A.5) 



Differentiating this equation with respect to z, dividing the result by d 2 A and applying 
d z leads to: 

dlA = d 2 z A[e d- z A + e 1 ], d 3 z A = d 2 z A[e 2 d z A + e 3 ], e , e u e 2 , e 3 e TZ . (A.7) 



Substituting these results into the equation (|A.5|) yields for A(z) and A(z) the differential 
equations: 

d\A = l -{d- z A) 2 ev + e,d- z A + e 5 , d 2 z A = ^(d z A) 2 e + e 4 d z A + e 5 , (A.8) 
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with the same constants eo,e4,es for both functions. Inserting these expressions into 
relation (|A.4j ) we get for dfipo = (m 2 +v)/2: 

e d 2 z AdlA 



2[f(d z A + d- z A) + e 4 } 2 



From equation ( |A.3| ) it follows that 
m 2 (l) = 

v(l) = 



!-?±)&AdlA 



^(d z A + d- z A) + e, 



1 + 



e Q 



2 J 



\ d 2 Adt 



A 



e 



(d z A + d- z A) + e 4 



1 -2 



(A.9) 

(A.10) 
(A.11) 



Thus we have to solve the differential equations (|A.7 ) for d z A and d z A with different 
choices of the constants e , e 4 and e 5 . Depending on the sign of A = 2e e 5 — e\ there exist 
three different cases: 



IV.l) A = 
IV.2) A < 



IV.3) A > 



d z A 
d z A 
d- z A 
d z A 



e 
A 



tan 



c(z - z 
/ A 



, d- z A 



A 



e 



tan 



2 
A 



V^A 



e 



<9 2 A = - 



v^A 

e 



tanh 



tanh 



e4 

£4 

(* - z ) I - — , 
J e 



c(z - Z ) 



(z - Zq) - — 



V^A 



2 
A 



(z - Zq) ) - j- 



J e 



For the case IV.l) we obtain: 



m 2 {l) 



4 



1 - 



e 



P-/o]- 



1 + 



eo 

2 ; 



20,^0 e^A.12) 
(A.13) 
(A.14) 
(A.15) 
(A.16) 



l = z + z . (A.17) 



If we choose eo = 2, then m (I) vanishes and v = — 2[l — la] is for Iq = identical with 
9^y(y?o) i n Liouville's model and on the other hand if one would like to obtain the 
PDE (|6.3|) only (this means m 2 (l) = — 2[l — /o]~ 2 ) one has to calculate eo from the eq. 
—2 = 4(1 — e /2)/eQ yielding e 0l2 = —1,2. Therefore one obtains the wave equation or 
once again Liouville's model. The last case represents only a map of (|6.3| ) onto itself (for 
details of Auto-BTs see e.g. |§). 

The case IV.2) yields: 



m 2 (l) 



A 

4 



i - 



e 



sin 



A 



1 + 



e 



sin 



{I ~ /o) 



(A.18) 

This is of no interest for our special models. Notice, however, that the special choice 
v = — A2 _1 sin _2 (vA[Z — lo]/2) leads to the wave equation. 

The case IV.3) is obviously similar to IV.2). The resulting functions m 2 (l) and v(l) are: 



A 



( o 



e 



m (I) — —z ( 1 — — ) sinh 



V^A 



il-k 



A 



( o 



-r 1 + — sinh 



v^A 



(A.19) 



38 



As the equations ( |6.11 ) and (|6.12 ) for the 4 models are not contained in the cases 
discussed up to now we return to case I: d z A = ai and d z A = a 1 . Assuming m 2 to be 
independent of I equation ( |A.3| ) then gives: 



dfi[) + {di^of + («i + ai)diip + ai<2i - m 2 



0. 



(A.20) 



Substituting £ = <9;-?/>o + ( a i + a i)/2 and K = m 2 + {a\ — ai) 2 /4 leads to: 



1.1) A > : i = A tanh(A(7 - l )) , 

1.2) A < : £ = - Atan(A(/ - /„)) 

1.3) A = : £ ' 



u = 2A 2 (1 - tanh 2 [A(Z - l )}) - m 2 (A.21) 
v = - 2A 2 (l+tan 2 [A(/ - l ))) - m 2 A.22) 

' 2 



m 



(A.23) 



The case 1.1 is the essential one for us. Because of the special type of the equations 
( p . 1 1| ) and ( |6.12|) , we are able to eliminate A by a transformation of variables z — > Az 
and z Az. Thus only the class of PDEs with a coefficient 2 in front of their tanh (7): 



d z d s 6 = jn tanh 2 (/) + n | 8 , n = 2 , n E 1Z 



(A.24) 



can be reduced by one BT to a Klein-Gordon or wave equation, like the relations we 
obtain in the Sine-Gordon ( |6.5|) and the Sinh-Gordon (|6.6|) theory. For the </> 4 models 
with n = 6 we need at least two BTs. 

Substituting m 2 ( y l) = 2ip (l)—v(l) in equation (|A.3|) leads to: 

d?*Po - Wo + ai)($Vo + SO - u = . (A.25) 

Choosing u to be equal to —brf{l) — c where the function r\ has to obey the equation 
dir) = brj 2 + c with b,c ElZ, making the special ansatz diipo = do + dir) with two constants 
do,di, inserting all this into equation ( |A.25| ) and comparing the coefficients in front of 
the powers of rj yields: 



rj 2 : d\b + b — d\ = , d\ 12 
rf_ : 2d + a\ + d\ = 

7?° 



± 



2c?o — — Oi — Oi , 

• ai = ±.2\l —die — c + ai. 



(A.26) 

(A.27) 
(A.28) 



dic + c - (d + a 1 )(d + a x ) = , 

For m 2 (Z) we obtain: 

m 2 {l) = 2dftp - v = c + 2d x c + {2d x b + b)rf . (A.29) 
Of special interest is the coefficient in front of rj 2 . Inserting d\ into eq. ( |A.29| ) we get: 

(A.30) 



2d x b + b = b 2 + b± b\/b 2 + 4b, 
with b = n(n + l) and choosing b to be equal to ±1 we obtain: 



2dtb + b = n{n + 1) + 1 ± ^1 + 4n(n + 1) 
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n(n — 1) 

(n + l)(n + 2) 



(A.31) 



Thus starting with a coefficient n(n + l) one BT can raise or lower n by 1. Choosing 
n = 1, 2, . . . we are able to calculate the functions ipi, which are essential in order to 
determine the i th BT of the hierarchy of n BTs, using the equation ( |A.2| ). Moreover 
the solutions ipi = exp(— ipj) of the i th PDE d z d z 6 + Vi$ = which reciprocal fulfils 
d z d z 9 + mj(l)9 = (see (|5.9|)) can be calculated. Obviously the special case b = —l and 
c = 1 leads to: 

d z d- z 9 = [n(n + 1) tanh 2 (^ + z) + a)B . (A.32) 

For n = 2 and special choices of a this relation yields the PDEs ( |6.11|) , (|6.12|) for the 
4 -theories. So they can be solved by n = 2 BTs. 
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Figures 



Here we display the 2-dimensional projections of the 2-dimensional extremals ( Q7.45 )) 
on the plane l = ^(t — vx) = const, and those 1-dimensional wave fronts ( |8.41| ) that lie 
entirely in this plane. The "singular" case (|8.41|) is shown in Fig. 1, the "regular" ones 
for c = —2 and c = 1 in Figs. 2 and 3, respectively. The extremals are plotted 
in solid lines, whereas the wave fronts are given in dotted ones. The axes are the field 
variable and the independent variable l=j(x — vt) which parametrizes the kink solution, 
obtained from the static one by a Lorentz transformation (see (|6.9|)). 



Fig. 1: The singular case 




In this case the transversality relations are fulfilled in the extended configuration space 
M.2+1 except on the planes (p = 2kn, k = 0, ±1, ±2, . . ., where the extremals and wave 
fronts are parallel. 
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Fig. 2: The regular case for c = — 2 




Here the transversality relations are fulfilled outside the planes <p = 2kir ± tt/3, k = 
0, ±1,±2, where the extremals and wave fronts are parallel since the Lagrangian 
density vanishes there. Notice, that the transversality relations are fulfilled on the 
boundaries of the regions M.^ = \2k-K < ip < 2{k + z, z G 1Z}, where the families 
of extremals ( f7.45|) are defined. Therefore the wave fronts can be continued from one of 
these regions to the next smoothly, contrary to the extremals. 
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Fig. 3: The regular case for c = l 



12 






10- 
8- 


****•• / / / / 

/ / 7 7' 


/////■ '/■ '// 

'*/*//// //**** 


6 


*' 3 T 4 ''//'/ 




4- 






2 






o- 






1 , , , , 1 





Notably, not only the transversality relation is satisfied everywhere in the extended 
configuration space, but in addition the wave fronts intersect the extremals everywhere 
perpendicularly. 
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